arXiv:1501.01101v2 [astro-ph.EP] 7 Jan 2015 


To APPEAR IN The Astrophysical Journal 
Preprint typeset using IATeX style emulateapj v. 5/2/11 


VALIDATION OF TWELVE SMALL KEPLER TRANSITING PLANETS IN THE HABITABLE ZONE 

Guillermo Torres^’^, David M. Kipping^ Francois Fressin\ Douglas A. Caldwell^, Joseph D. Twicken^, Sarah 
Ballard^’^, Natalie M. Batalha®, Stephen T. Bryson®, David R. CiardT, Christopher E. Henze®, Steve B. 
Howell®, Howard T. Isaacson®, Jon M. Jenkins®, Philip S. Muirhead^’^®, Elisabeth R. Newton^ Erik A. 

Petigura®, Thomas Barclay®, William J. Borucki®, Justin R. Crepp^\ Mark E. Everett^^, Elliott P. Horch^®, 
Andrew W. Howard^^, Rea Kolbl®, Geoffrey W. Marcy®, Sean McCauliff^®, and Elisa V. Quintana® 

To appear in The Astrophysieal Journal 

ABSTRACT 

We present an investigation of twelve candidate transiting planets from Kepler with orbital periods 
ranging from 34 to 207 days, selected from initial indications that they are small and potentially in 
the habitable zone (HZ) of their parent stars. Eew of these objects are known. The expected Doppler 
signals are too small to confirm them by demonstrating that their masses are in the planetary regime. 

Here we verify their planetary nature by validating them statistically using the BLENDER technique, 
which simulates large numbers of false positives and compares the resulting light curves with the 
Kepler photometry. This analysis was supplemented with new follow-up observations (high-resolution 
optical and near-infrared spectroscopy, adaptive optics imaging, and speckle interferometry), as well as 
an analysis of the flux centroids. Eor eleven of them (KOI-0571.05, 1422.04, 1422.05, 2529.02, 3255.01, 

3284.01, 4005.01, 4087.01, 4622.01, 4742.01, and 4745.01) we show that the likelihood they are true 
planets is far greater than that of a false positive, to a confidence level of 99.73% (3cr) or higher. 

Eor KOI-4427.01 the confidence level is about 99.2% (2.6cr). With our accurate characterization of 
the GKM host stars, the derived planetary radii range from 1.1 to 2.77^0. All twelve objects are 
confirmed to be in the HZ, and nine are small enough to be rocky. Excluding three of them that have 
been previously validated by others, our study doubles the number of known rocky planets in the HZ. 
KOI-3284.01 (Kepler-438 b) and KOI-4742.01 (Kepler-442 b) are the planets most similar to the Earth 
discovered to date when considering their size and incident flux jointly. 

Subjeet headings: methods: statistical — planetary systems — stars: individual (KOI-0571 = Kepler- 
186, KOI-1422 = Kepler-296, KOI-2529 = Kepler-436, KOI-3255 = Kepler-437, 
KOI-3284 = Kepler-438, KOI-4005 = Kepler-439, KOI-4087 = Kepler-440, KOI- 
4427, KOI-4622 = Kepler-441, KOI-4742 = Kepler-442, KOI-4745 = Kepler-443) 

— techniques: photometric 


1. INTRODUCTION 

Over the duration of its four-year mission the Kepler 
spacecraft has enabled the identification of several thou¬ 
sand candidate transiting planets {Kepler Objects of In¬ 
terest. or KOIs: [Borucki et al.ll2Qlll: iBatalha et al.ll2QT^ 
iBurke et al.ll2QlJ ). and many more continue to be found 
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from reanalysis of the original data with increasingly so¬ 
phisticated methods. Only a tiny fraction of these can¬ 
didates have been “confirmed” in the traditional sense 
of having had their masses measured, either spectro¬ 
scopically or by modeling their transit timing variations 
(TTVs). Hundreds of others, mostly in multiple systems 
(“multis”), have recently been shown stati stically to have 
a very high chance of b eing true planets (jLissauer et al.l 
12(111 IBnwe ^nill2(rf4l . even if their masses are not cur¬ 
rently known. 

Public curiosity and scientific interest have motivated 
efforts in the last few years to find and confirm rocky 
planets similar in size to the Earth that are orbiting in 
the so-called habitable zone (HZ) of their parent stars, 
usually taken in this context to be the region in which 
water on t he surface can be in a liquid state. Recent es¬ 
timates by iPetignra et al.l (j2Q13[ ) suggest the rate of oc¬ 
currence of Earth-size planets in the HZ of Sun-like stars 
(which they defined broadly as corresponding to an inci¬ 
dent flux between | and 4 times that of the Earth) may 
be as high as 22 T 8 % (fo r other estimates seelZsom et al.l 
l2Q13l : iPintr et ahl l2Ql4 lEoreman-Mackev et al.l l2Q14f ). 
but very few of these have actually been confirmed. 
These two conditions for habitability — rocky nature 
and suitable location — are sometimes more challenging 
to establish unambiguously than it may seem, as they 
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require an accurate knowledge of the host star prop¬ 
erties, something that is not always trivial to achieve 
for the faint Kepler targets. Definitions of the HZ 
are actively debate d and have changed over time (e.g 


Kasting et al.lll99^ iSelsis et al.l 200 

iLeconte et al.ll 2 Ql 


20131: ISeagerll20l4 ID 

Kopparapu et al.H2Q14f ). as has our understanding of the 
radius at which planets transition from being rocky to 
being dominated by thick hydrogen/helium envelopes, 
which presumably makes them unsuitable for life as 
we know it. The interface i s currently thought to 


Kopparapu et ail 

Zsom et al.ll2Q13l: 


be between 1.5 and 2 Ra, fe.g.. iLopez fc 


Weiss fc Marcvl [2Ql4 llV^rcv et al.l l2Ql4 


ortne 


Rogeri 


i 


2013; 


mm 


Wolfgang fc Lopezll2014f ) , but many uncertainties remain 


and it is possible that some planets that are slightly 
larger may still be rocky. 

These difficulties notwithstanding, a handful of tran¬ 
siting planets that appear to meet these conditions have 
already emerged from the Kepler sample. Among those 
considered at the time of publication to be smaller than 
2.5 i ?0 and in or near the HZ, the firs t to be announced 
was Kepler-22 b (jBorucki et al.l[ 20 T^ . with a radius of 
Rp = 2.38 ± O.13i^0. Others followed incl uding Kepler 


61b (2.15 ± O.13i^0; iBallard et all [20T^. Kepler-62e 
and f (1.61 ± 0.05 and 1.41 ± 0.07 R^; Borucki et al.l 
l2013[ ). Kepler-69c (I.TI^q oq i? 0 : l^arclav et~al] l2013l ). 
and r ecently Kepler-186 f (1.11 TO. 14770 : [Quintana et al.l 
l2014i ) . All of these were ‘Vali dated” statistically using 
techniques such as BLENDER ([Torres et al.l 120041 1201 ll 
iFressin et al.l [20T^ by examining the likelihood of dif¬ 
ferent false positive scenarios to measure the confidence 
in the planetary interpretation. Some of them are in 
multi- planet systems; as mentioned above. iLissauer et al.l 
(j2014f ) have developed a statistical argument by which 
they demonstrate that most candidates in multis are 
very likely to be true planets. Based on this reason- 
ing and a care f ul exa mination of follow-up observations, 
ILissauer et al.l (j2014f ) reported an additional small planet 
(Kepler-296 f, 2.3 1 T 0.12 77^) that appears to be in the 
HZ, according to iRpwe et al.l (j2014f ). Taking advantage 
of the same statistical framework for validation, the lat¬ 
ter authors announced one more small HZ planet in 
a multi-planet system: Kepler-174d (2.19 ± 0.13 770 ). 
Other small HZ candid a tes have been proposed (see, e.g., 
Muir head et al.l 120121: [Dressing <V Charbonneaul 2013 : 


GaidosI 1^^ Mann et al.l l2013bl: iPetigura et al.l 2013 : 
Star et al.ll2014l) but have vet to be investigated in detail 


and confirmed. A few additional examples also believed 
to be in this special group were later shown to fail one or 
both requirements (size and location), either because of 
incorrect stellar parameters (originally based on photom¬ 
etry and later improved with spectroscopy) or because of 
the discovery of other stars in the photometric aperture 
that dilute t he transits and c hange the inferred planetary 
radius fe.g.. lStar et al.l[2QT3 ). 

In this paper we investigate a sample of twelve Kepler 
candidates identified on the basis of preliminary analy¬ 
ses indicating that they are small and likely to be in or 
near the HZ of their parent stars. Our goal is to val¬ 
idate them as bona-fide planets, and to confirm their 
key properties pertaining to habitability. We describe 
our efforts over the past year to obtain the follow-up 
observations necessary to robustly characterize the host 


stars and to validate the signals. In the interim, three of 
them that are in systems with five candidates each (KOI- 
0571.05 = Kepler-186 f, KOI-1422.04 = Kepler-296 f, and 
KOI-1422.05 = Kepler-296e) were the subject of re¬ 
cent studies by others. KOI- 1422.04 and KOI-14 22.05 
were annou n ced a s planets by ILissauer et al.l (j2014l ) and 
iRowe et al.l (j2014f ) using their statistical framework for 
validation of multis “in bulk”, though only the first was 
mentioned as bein g in the HZ. KO I -0571 .05 was validated 
independently by lOuintana et al.l (j2014l) . These authors 
relied in part on the fact that false positives are much less 
common in multis, and on the validation of the other four 
planets in the same system by the previous authors, in or¬ 
der to reach a sufficiently high level of confidence for the 
fifth candidate. Most of the other targets in our sample 
are considerably more difficult to validate because they 
tend to have long orbital periods and do not belong to 
multiple systems. Nevertheless, as we show below, the 
application of the BLENDER technique used in many of the 
previous discoveries is able to achieve validations here as 
well. Our work has now essentially doubled the number 
of known small HZ planets. 

Our paper is organized as follows: Sect. [ 2 ] and Sect. [3] 
describe our sample of candidates and the Kepler pho¬ 
tometry we use. The follow-up observations including 
high-resolution imaging, optical spectroscopy, and near- 
infrared spectroscopy are presented in Sect. IH where we 
also report the complicating discovery of nearby stellar 
companions to four of the targets. These close compan¬ 
ions not only dilute the transits but also bring ambiguity 
as to the precise location of the planets in these sys¬ 
tems. A description of the analysis of the ffux centroids 
that addresses some of those companions is given here as 
well. Sect. [5] follows with a determination of the stellar 
properties. The formal validation of our candidates is 
the subject of Sect. [ 6 l after which we proceed with the 
transit light curve fits to derive the planetary parameters 
(Sect. [7j). In Sect. [ 8 ] we apply the powerful technique of 
asterodensity profiling to investigate whether the planets 
that are in multiple systems orbit the same star, to ex¬ 
tract eccentricity information from the light curves of our 
targets, and for the four host stars with close stellar com¬ 
panions, to address the possibility that the planets orbit 
the companions rather than the intended targets. The 
topic of habitability is examined in Sect. [9l We conclude 
with a discussion of our findings and final remarks. 

2. CANDIDATE IDENTIFICATION 

Candidate transiting planets of special interest for 
this study were initially identified as ‘threshold crossing 
events’ (TCEs; objects with a 7.1a or greater multiple- 
event transit detection statistic; see I Jenkins et al.ll 2 QQ 2 L 
| 2 Q 1 Q[) from a transit search conducted using nearly three 
years of Kepler data gathered from quarters 1 to 12 , with 
version 8.3 of the Science Operations Cen ter (SOC) pro¬ 
cessing pipeline (jTenenbaum et aP l 2 nn . The analysis 
of three years of data gave us the first chance to de¬ 
tect three transits from planets in habitable-zone orbits 
around stars like the Sun, as well as to significantly in¬ 
crease the signal-to-noise ratio (SNR) for transit signals 
from habitable-zone orbits around cooler stars. In order 
to investigate the smaller planets potentially in the HZ, 
we selected for further scrutiny the TCEs that met the 
criteria Rp < 2 R^ and Teq < 303 K, as well as a transit 
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SNR > 7. Here Teq represents the planetary equilibrium 
temperature, calculated assuming an Earth-like albedo, 
and the stellar properties r elied on informatio n from the 
Kepler Input Catalog (KIC: [BFown et al.1l2Qll[) . We used 
data products from versions 8.3 and 9.0 of the SOC Data 
Validation (DV) pipeline module for the initial triage and 
vetting (for an overview of this process, see IWu et ahl 
[2oT(I FI. We compared the results of our vetting with the 
scores from early runs of a n automated vetting c ode be¬ 
ing developed at the time (jMcCanliff et al.ll2Q14f ). From 
the 18,407 TCEs that passed through DV, we identified 
385 TCEs around 338 unique stars that met our selec¬ 
tion criteria. From this list of 385 TCEs we identified ten 
that passed fiux and centroid vetting criteria and that 
had low false-alarm probabilities from the auto-vetter. 
During this vetting process, the SOC completed a tran¬ 
sit search of Kepler data from quarters 1-16 using ver¬ 
sion 9.1 of the processing pipeline. With the addition of 
three quarters of data and the improvements in the ver¬ 
sion 9.1 DV products, we promoted an additional TCE 
(KOI-4742.01), which had marginal signal-to-noise ratio 
in the 12-quarter run, and identified KOI-1422.05, the 
fifth KOI and second habitable-zone candidate around 
this star. The complete list of targets selected for this 
study is given in Table [TJ 

We note that there is significant contamination of the 
TCE list for periods near the orbital period of the Ke¬ 
pler spacecraft (372 days) caused by instrumental im¬ 
age artif acts on specific detector channels, as has been 
noted by iTenenbanm et all (j2Q13f ) . The contamination 
is due to moire pattern noise inj ected during readout 
of some of the detector chains (jCaldwell et al.l I2Q1QI : 
iKolodzieiczak et ^ I2Q1Q[ ) . The moire noise generates 
a host of single-event transit-like signals for targets on 
the noisy channel. When folded at a period near one Ke¬ 
pler year these signals can give multiple-event statistics 
above the 7.1a detection threshold, as the same target 
star falls on the same noisy channel. We have not in¬ 
cluded any KOIs from the moire pattern noise channels 
in our analysis here. Two of the KOIs (KOI-2529.02 and 
KOI-4005.01) do show some of their transits on detec¬ 
tor channel 56, which has increased noise relative to the 
other channels on which these targets fall. However, we 
included these KOIs in our sample because both show 
transits on other channels and both passed all of our 
vetting checks. 

3. KEPLER PHOTOMETRY 

For the analysis in this paper we have made use of 
the publicly available Kepler data for each KOI from the 
Mikulski Archive for Space Telescopes (MAST)E3 The 
observations span Aep/er quarters 1-17, corresponding to 
a period of four years beginning in 2009 May. The down¬ 
loaded data were made available as part of Data Release 
23 and were processed using SOC Pipeline version 9.1. 
While long-cadence (29.4 min) data were available for all 
KOIs, short-cadence (58.9 s) data were only available for 
KOI-0571 and KOI-1422. In the following we describe 
our further processing of the simple aperture photome¬ 
try (SAP) time series used for the detailed light-curve 

All DV results are available from the NASA Exoplanet 
Archive: http: //exoplanetarchive . ipac. caltech. edu . 

https://archive.stsci.edu/index.html . 


analysis of all KOIs. Later in Sect. 15.11 and Sect. 16.11 we 
describe the slightly different processing of the photom¬ 
etry used to extract stellar rotation information and for 
the BLENDER analysis. 

Prior to making use of the SAP measurements we re¬ 
moved instrumental and stellar photometric variability 
in the Kepler data that can disturb the transit light 
curve profile. This process was conducted in two steps: 
pre-detrending cleaning, and long-term detrending. The 
pre-detrending cleaning was carried out independently 
for each transit event, restricted to plus or minus half 
an orbital period surrounding the mid-transit time. We 
visually inspected each epoch and removed any charge¬ 
trapping ramps, flare-like events, and instrumental dis¬ 
continuities in the data. We made no attempt to cor¬ 
rect these artifacts and simply excluded them from the 
photometry manually. We then removed all of the tran¬ 
sit signals of the other known candidates in each system 
(for stars having multiple candidates) within ±0.6 transit 
durations centered on the mid-times, using the reported 
ephemerides. Finally, we cleaned the data of 3a outliers 
from a moving median smoothing curve with a 20-point 
window. 

Next we removed long-term trends, which can be due 
to instrumental effects such as focus drift, or stellar ef¬ 
fects such as rotational modulations. For this task we 
used the Cosine Filtering with Autocorrelation Mini¬ 
mization (CoFiAM) algorithm, which was specifically de¬ 
veloped to protect the shape of a transit light curve. 
We briefly highlight the main features of CoFiAM and 
the specific input s used for this analy sis, and we di¬ 
rect the reader to iKipping et al.l (j2Q13f ) for a more de¬ 
tailed description of the technique. It is essentially a 
Fourier-based method that removes periodicities occur¬ 
ring at timescales greater than a pre-designated “pro¬ 
tected timescale”. In this work, we selected three times 
the known transit duration to be the protected timescale, 
which ensures that the transit profile is not distorted 
in frequency space. CoFiAM does not directly attempt 
to remove high frequency noise, since the Fourier trans¬ 
form of a trapezoidal-li ke light curve con tains significant 
high frequency power (jWaldmannI l2Q12[ ). It is able to 
explore many different harmonics by trying longer pro¬ 
tected timescales than the nominal choice (we capped 
the maximum number of harmonics at 30) and evalu¬ 
ate the autocorrelation at a pre-selected timescale (we 
used 30 minutes) locally surrounding each transit. From 
these different harmonics, we selected the harmonic order 
that minimizes this local autocorrelation, as quantified 
using the Durbin-Watson statistic. This “Autocorrela¬ 
tion Minimization” component of CoFiAM provides opti¬ 
mized data for subsequent analysis. For each KOI we 
defined the local transit data as being within six tran¬ 
sit durations either side of the mid-transit time, in order 
to provide an ample out-of-transit baseline. These local 
data were divided through by the final CoFiAM function 
and then stitched together to form our final light curve 
for analysis. 

4. FOLLOW-UP OBSERVATIONS AND CENTROID 
ANALYSIS 

4.1. High-resolution imaging 

The presence of other stars blended with our targets 
is a potential concern, as these companions could be or- 
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Table 1 

Sample of KOIs in this study. 


Candidate 

Name 

KID 

Kp 

(mag) 

b 

(deg) 

Period 

(days) 

Depth 

(ppm) 

KOI-0571.05 

Kepler-186 f 

8120608 

14.625 

8.2 

129.94 

540 

KOI-1422.04 

Kepler-296 f 

11497958 

15.921 

18.0 

63.34 

1070 

KOI-1422.05 

Kepler-296 e 

11497958 

15.921 

18.0 

34.14 

850 

KOI-2529.02 

Kepler-436 b 

8463346 

15.856 

6.5 

64.00 

1150 

KOI-3255.01 

Kepler-437 b 

8183288 

14.352 

9.0 

66.65 

520 

KOI-3284.01 

Kepler-438 b 

6497146 

14.467 

18.6 

35.23 

400 

KOI-4005.01 

Kepler-439 b 

8142787 

14.560 

19.9 

178.14 

700 

KOI-4087.01 

Kepler-440 b 

6106282 

15.134 

15.8 

101.11 

1010 

KOI-4427.01 


4172805 

15.645 

8.2 

147.66 

1220 

KOI-4622.01 

Kepler-441 b 

11284772 

15.142 

19.1 

207.25 

960 

KOI-4742.01 

Kepler-442 b 

4138008 

14.976 

15.0 

112.31 

560 

KOI-4745.01 

Kepler-443 b 

11757451 

15.891 

17.0 

177.67 

930 


Note. — Columns after the first indicate the Kepler planet designation, 
Kepler identification number, brightness in the Kepler passband. Galactic 
latitude, orbital period, and transit depth in parts per million relative to the 
out-of-transit stellar fiux. For consistency in this paper we will refer to all 
objects by their original KOI names throughout, even though the validations 
described later earn them the official Kepler planet designations listed in the 
second column. 


bited by another object causing eclipses that may be the 
source of the transit signals we observe. Even if they are 
not eclipsing, the light of the companion stars will atten¬ 
uate the signal and lead to a bias in the planetary radius 
determination if the effect is not accounted for. 

Images fro m the J-band UK In frared Telescope sur¬ 
vey (UKIRT dLawrence et al.l[2QQ7[) have indeed revealed 
companions to three of our targets (KOI-2529, 3284, and 
4427) that are close enough to fall within the photo¬ 
metric aperture of Kepler (i.e., within a few arc sec¬ 
onds), although as we describe later they are ruled out 
as th e so urce of the transits by our centroid analysis in 
Sect. 14.21 We list these companions in Table [21 based on 
information taken from the Kepler Community Follow¬ 
up Observing Program (CFOP) Web sit43. Additional 
companions to three other KOIs are seen in the UKIRT 
images at larger separations (and are also excluded by 
the centroid analysis), but these stars are all very faint 
(Am > 5), and have a negligible impact on the inferred 
planetary sizes. 

The UKIRT images have a typical seeing-limited res¬ 
olution of about 078 or 079. To explore the inner re¬ 
gions around our targets beyond the reach of UKIRT we 
observed them with near-infrared adaptive optics (AO) 
in the J (1.248 um) and K' (2.124 jum) filters using the 
NIRC 2 imager (jWizinowich et al]l2004 [Johansson et ah! 
[200l on the Keck II, 10 m telescope. KOI-0571 and 
KOI-3255 were observed in August 2012, and the rest 
in August 2013 a s part of a general infrared AO sur¬ 
vey of KOIs (e.g. , [Adams et al.[ \2012l [Rowe et all [2014 
[Marcv et al.[[20l4 . KOI-1422 was not observed with AO 
for this project a s it has been the target of a separate ef - 
fort using HST (|Gilliland et al.[[2015 : [Star et al.[[2014D : 
those authors found it to have a close companion (see 
below). For KOI-0571, 3255, 2529, 3284, and 4742 the 
targets themselves were used as natural guide stars; for 
KOI-4005, 4087, 4427, 4622, and 4745 we used the laser 
guide-star AO system. In all cases the observations were 
obtained in a 3-point dither pattern to avoid the lower 
left quadrant of the NIRC2 array, which displays elevated 

https://cfop.ipac.caltech.edu/home/ . 


Table 2 

Close companions to target stars. 


Star 

Aug. sep. 
(") 

P.A. 

(deg) 

Mag. diff.^ 
(mag) 

Source 

KOI-1422 

0.220 

216.3 

AR 

= 1.72b 

Speckle 

KOI-2529 

0.402 

350.4 

AK‘ 

' = 4.69 

Keck AO 


5.05 

42.6 

AJ-- 

= 3.75 

UKIRT 

KOI-3255 

0.180^ 

336.5 

AK 

' = 0.05 

Keck AO 

KOI-3284 

0.438^ 

193.2 

AK 

' = 2.03 

Keck AO 


4.10"^ 

359.7 

AJ - 

= 2.90 

UKIRT 

KOI-4427 

4.76 

274.8 

AJ : 

= 2.87 

UKIRT 


^ For targets with both a close and a wide companion the 
magnitude differences for the wide companions are relative 
to the total brightness of the inner pair. 

^ Also observed at 880nm (approximately Sloan 2 : band), 
giving A 2 : = 1.62. The separation and P.A. in the table are 
the average of the two s peckle bands. Similar results were 
reported bv IStar et al.l 1201411 from HST observations. 

^ This companion was also detected in our speckle imaging 
observations. 

^ This companion was also detected in the UBV survey 
of lEverett et al U 20 ID . with magnitude differences AB = 

1.80 and AV = 2.01. The separation and P.A. in the table 
are the average of the two surveys. 

noise. Five images were collected per dither pattern po¬ 
sition, each shifted 075 from the previous dither position 
to enable the use of the source frames for the creation of 
the sky image. 

In general the NIRC2 imaging detects all of the sources 
found in the UKIRT J-band imaging within 5" of each 
target. The NIRC2 array has 1024 x 1024 pixels with 
a scale of about lOmaspix”^, and a field of view of 
1071 X 1071. Each frame was dark-subtracted and fiat- 
fielded, and the sky frames were constructed for each 
target from the target frames themselves by median fil¬ 
tering and co-adding the 15 or 25 dithered frames. In¬ 
dividual exposure times varied depending on the bright¬ 
ness of the target, but were long enough to obtain at 
least 5000 counts per frame (NIRC2 has a gain of 4 elec¬ 
trons per DN); frame times were typically 1-30 seconds. 
Data reduction was performed with a custom set of IDL 
routines. Close companions were detected around KOI- 
2529, 3255, and 3284 (see Figure [T]), two of which also 

































Validation of small HZ planets 


5 



1.00 

0.50 

0.0 

-0.50 

- 1.00 

1.00 0.50 0.0 - 0.50 - 1.00 

Aa [axcsec] 



1.00 

0.50 

1 
*0 
<1 

-0.50 

- 1.00 

1.00 0.50 0.0 - 0.50 - 1.00 

A<x [arcsec] 



Figure 1. Images in the K' band from our Keck AO imaging 
observations showing close companions found around KOI-2529, 
KOI-3255, and KOI-3284 (see also Table [2|). 



Figure 2. Sensitivity curves from our Keck AO imaging observa¬ 
tions in the (2.124/rm) K' band. Curves correspond to the KOIs 
as labeled, top to bottom. 


have wider companions seen in the UKIRT images. We 
report the relative positions and brightness of these close 
companions also in Table [2l 
Point source detection limits were estimated in a se¬ 
ries of concentric annuli drawn around each target star. 
The separation and widths of the annuli were set to the 
full width at half maximum (FWHM) of the target point 
spread function (PSF). The standard deviation of the 
background counts was calculated for each annulus, and 
the 5cr li mits were determin ed for each annular region 
(see also lAdams et al.1l2Ql2f ). The PSF widths for our 
Keck images were typically found to be 4-5 pixels, cor¬ 
responding to 0'.'04-0'.'05 FWHM. Typical contrast levels 
from the images are 2-3 magnitudes at a separation of 1 


FWHM, 7-8 magnitudes at 5+ FWHM, and deeper past 
10 FWHM. Sensitivity curves in the K' band for each 
of the targets are shown in Figure [21 and extend to an¬ 
gular separations of approximately 3", set by the dither 
pattern overlap. We supplemented these measurements 
with similar sensitivity estimates made from the UKIRT 
images available on the CFOP Web site, which reach to 
much wider separations but are not as deep. 

Additional high-resolution imaging observations for 
four of our targets (KOI-0571, 1422, 3255, and 3284) 
were obtained wi th the Differential Speckle Survey In- 
strument (DSST: iHorch et ^l2QQ9l : (Howell et al.ll2Qlll ) 
mounted on the 8.1m Gemini-N telescope during the 
time period UT 2012 July 25-28. The seeing conditions 
were generally 075-078 during those nights. The large 
aperture of the Gemini-N primary mirror allows for high 
SNR observations at faint magnitudes and achieves a 
very small diffraction limit for visible-light observations. 
The speckle camera provides simultaneous observations 
in two filters. For all Gemini-N observations discussed 
here we used filters centered at 692 nm and 880 nm (cor¬ 
responding very roughly to the Gousins R and Sloan 2 ; 
bands), with band widths of 40nm and 50nm, respec¬ 
tively. Although full details of the observing prot ocols 
have been given previously by iHorch et ah! (|2012f ) , we 
summarize them here for completeness. 

For bright targets observed with this instrument (U < 
12) a single sequence of 1000 speckle frames of 60 ms du¬ 
ration each has usually been found to be sufficient, for a 
total of 3-4 minutes of observing time. The much fainter 
sources discussed here (15 < V < 17) required up to 
eight, 1000-frame sequences that were later combined, or 
about 30 minutes of on-source exposure time each. Cal¬ 
ibration images and sequences were taken as is usual to 
allow us to measure and define dispersion effects, position 
angles, faint limit sensitivity, photometric performance, 
and for point source reconstruction. 

The speckle observations revealed no stars in the vicin¬ 
ity of KOI-0571, but KOI-1422 was found to have a 
close companion (seen in both filters) that was also 
subsequently detecte d in the RST observations by 
I Gilliland et ah! (j2Q15[ ). Its brightness and position rel¬ 
ative to the primary are given in Table [21 Close compan¬ 
ions were also found around KOI-3255 and KOI-3284, 
which are the same as detected in our AO imaging with 
Keck. 

For estimating our detection threshold for faint com¬ 
panions we considered each local maximum in the re¬ 
constructed image as a potential stellar source, and de¬ 
termined the statistics of these peaks (i.e., their average 
values and standard deviation) as a function of angu¬ 
lar distance from the primary star. We then adopted a 
conservative 5a threshold for the detection of any com¬ 
panion stars. The details of the calibration procedures 
mentioned abov e and the precision that rn ay be obtained 
are described by IHorch et al.l (j2QllL l2Q12f ) . The sensitiv¬ 
ity curves for the four Kepler targets observed here with 
speckle imaging are shown in Figure O 

4.2. Centroid analysis 

One method to identify possible false positives due to 
background eclipsing binaries measures the location of 
the transit signal relative to each KOI host star via dif¬ 
ference imaging. For each quarter, a difference image 
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Figure 3. Similar to Figure [2] for the four targets observed with 
speckle interferometry at the Gemini North telescope. Curves are 
shown for the (692 nm) R-hand observations, and correspond to 
the KOIs as labeled, top to bottom. 

may be formed by subtracting an average of in-transit 
pixel values from out-of-transit pixel values. If the tran¬ 
sit signal is due to a stellar source, the difference image 
will show that stellar source, whose location is deter- 
mined by Pixel Res ponse Function centroiding (PRF; see 
iBrvson et al.ll2QlQ[) . This location may then be compared 
to the location of the host star, determined from a PRF- 
fit centroid of an average out-of-transit image (assuming 
the host is well isolated). Repeating this for each quar¬ 
ter in which transits occur gives a collection of offsets 
of the transit source from the host star. The average of 
these quarterly offsets defines the final offset of the tran¬ 
sit source from the host star. This metric is reported for 
each KOI on the NASA Exoplanet Archive. If the offset 
is greater than three times its uncertainty, there is cause 
for concern that the signal may be coming from a back¬ 
ground source rather than the KOI itself. In that case we 
would conclude that the likelihood that the transit signal 
is due to a planet around the KOI host star is low (for de¬ 
tails, see lBrvson et al.l[2Q13[ ). Any transit source location 
less than 3a from the host star is considered statistically 
indistinguishable from the host star. This 3a exclusion 
radius depends on the KOI rather than the host star, as 
it is primarily a function of the signal-to-noise ratio of 
the KOFs transit signal. 

Table [3] lists the average offset of the transit source for 
each KOI from the KOI host star, as reported on the Exo¬ 
planet Archive, as well as the 3cr exclusion radius for each 
KOI. All KOIs have offsets that are less than three times 
their uncertainties, with the exception of KOI-4427.01. 
Therefore, on centroid grounds alone KOI-4427.01 has a 
non-negligible chance of being a false positive. 

The companions at 4" to 5" from KOI-2529, 3284, and 
4427 listed in Table [2] are significantly farther than 3cr 
from the respective observed transit positions, so they 
are ruled out as possible sources of the transit, as are 
any wider companions detected in the UKIRT images. 
On the other hand, the close (< 075) companions to 
KOI-1422, 2529, 3255, and 3284 cannot be ruled out by 
centroid analysis alone. We address this issue later. Ad¬ 
ditional closer (unresolved) companions may of course 
also be present. 


Table 3 

Centroid results for the candidates. 


Candidate 

Offset from 
star (") 

Offset from 
star (cr) 

3f7 exclusion 
radius (") 

KOI-0571.05 

0.74 ± 0.82 

0.90 

2.45 

KOI-1422.04 

0.34 ± 0.34 

1.00 

1.03 

KOI-1422.05 

0.54 ±0.33 

1.67 

0.98 

KOI-2529.02 

0.08 ±0.64 

0.13 

1.92 

KOI-3255.01 

0.55 ±0.19^ 

2.82 

0.58^ 

KOI-3284.01 

0.60 ± 0.94^ 

0.64 

2.82^ 

KOI-4005.01 

0.23 ±0.61 

0.38 

1.84 

KOI-4087.01 

0.11 ±0.41 

0.26 

1.22 

KOI-4427.01 

0.83 ± 0.24^ 

3.47 

0.72^ 

KOI-4622.01 

0.92 ± 0.52 

1.76 

1.57 

KOI-4742.01 

1.19 ±1.17 

1.01 

3.52 

KOI-4745.01 

0.76 ±0.36 

2.09 

1.08 


^ The presence of relatively bright close companions to KOI- 
3255 and KOI-3284 can cause biases in the centroid measure¬ 
ments. These biases are expected to be smaller than the sep¬ 
aration of the close companion, but we have nevertheless con¬ 
servatively doubled the uncertainty found on the Exoplanet 
Archive. 

^ The Exoplanet Archive reports the offset of KOI-4427.01 as 
0794 ±1757, but there is a single large outlier due to a quarter 
with a noisy difference image. The values reported here are 
after removing the outlier. 

The centroid offsets in Table [3l on the basis of which 
we rule out the 4-5" companions, were computed on the 
assumption that the host stars do not have closer (< 4") 
companions of c omp arable brightness. However, as de¬ 
scribed in Sect. lO this is not the case for KOI-3255 
and KOI-3284. The relatively bright close companions 
for these two KOIs will introduce error in the centroid 
measurement of the KOI star position. (Such error does 
not arise for the difference-image centroids measuring the 
transit source position because, assuming negligible vari¬ 
ability, the difference image shows only the star hosting 
the transit signal.) To account for this unknown error, 
we have chosen to conservatively double the offset un¬ 
certainties for KOI-3255 and KOI-3284, although we do 
not expect the error to be larger than the separation of 
the companion. The uncertainties in the table already 
include this conservative doubling. The companion for 
KOI-2529 is sufficiently dim that the centroid error due 
to this companion should be negligible. 

4.3. Optical spectroscopy 

In order to characterize the target stars and search for 
additional companions, we acquired high-resolution spec¬ 
tra for the eleven stars in our sample with the Keck I 
Teles cope on Mauna K ea (HI) and the HIRES spectrom¬ 
eter (jVogt et al.l 119941) . using the s tandard setup of th e 
California Planet Search (CPS; see lHoward et al.l[2QlQ[ ). 
The wavelength coverage is approximately 3800-8000 A, 
recorded on three contiguous CCD detectors. Obser¬ 
vations were made over the course of four nights (UT 
2012 July 26 and September 22, and 2013 June 7 and 
28) under clear skies, with a typical seeing of ^1". Ex¬ 
posure times of about 15 minutes resulted in signal-to- 
noise ratios of 12-25 per pixel at a mean wavelength 
near the Mg Ib triplet (5150-5200 A). The exception was 
KOI-3284, which was observed for 45 minutes yielding 
a SNR of 45 per pixel. All spectra were taken using 
the C2 decker (0786 x 14"), giving a resolving power of 
R ^ 60,000. The 14" long slit allows for subtraction of 
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the sky background that is superimposed on the stellar 
spectrum. Sky subtraction was performed by measuring 
the sky background light on a wavelength-by-wavelength 
basis. 

For the reduction of the raw spectra we used the stan¬ 
dard CPS pipeline. Images from each of the three HIRES 
CCDs were independently reduced and subjected to sky 
subtraction, flat fielding, and cosmic ray removal. The 
pixel columns at each wavelength were then summed, 
providing photon counts as a function of wavelength for 
each pixel. Consistent wavelength solutions were ensured 
by aligning a carefully chosen set of Thorium-Argon 
emission lines onto the same pixels at the beginning of 
each night’s observations. 

The radial velocity (RV) of each star was measured 
using the A-band and B-band telluric line features as 
wavelength fiducials. We measured the placement of the 
stellar absorption lines relative to these features, and 
then referenced them to stars of known radial velocity 
(jChubak et al.ll2QT^ . The final radial velocity measure¬ 
ments in the frame of the solar system barycenter are 
accurate to ±0.1kms“^, and are given in Table |4] with 
other details of the observations. 

To aid in eliminating some of the false positives that 
might be causing the transit-like signals in our targets, 
we searched each of our spectra for evidence of a sec¬ 
ond set of lines from another star possibly falling on 
the spectrometer slit during our observations (i.e., at 
angular separations smaller than about 0743, or half 
of the slit width). We first cross-correlated each spec¬ 
trum against a large library of stellar spectra obtained 
with the same instrument and setup (see Sect. © in or¬ 
der to identify the best match. After flux normaliza¬ 
tion and placement of the target star and library stars 
onto a common wavelength solution, we subtracted this 
best match from the observed spectra, and the residu¬ 
als were subsequent ly cross-correlate d against the same 
set of library stars (jKolbl et al.ll2Ql^ . A resulting peak 
in the cross-correlation function would be evidence of a 
second star contaminating the spectrum. For our targets 
this method is sensitive to physically associated as well 
as unrelated (background/foreground) companion stars 
so long as the relative RV of the two objects is greater 
than about 10kms“^. For smaller velocity differences 
the technique loses sensitivity due to imperfect subtrac¬ 
tion of the primary star. 

We detected no companions to any of the 11 stars 
in our sample down to about 1% of the brightness of 
the primary (corresponding to a magnitude difference 
AKp = 5) for targets brighter than Kp =15, and down 
to 2% {AKp ^ 4.2) for our fainter targets. This indicates 
that at least two of the close companions found in our 
Keck AO imaging that are bright enough to have been 
detected spectroscopically, those around KOI-3255 and 
KOI-3284, must have radial velocities within 10kms“^ 
of the main star, which strongly suggests they are phys¬ 
ically associated with the target. 

4.4. Near-infrared speetroseopy 

For two of our cooler targets (KOI-4087 and KOI-4622) 
we gathere d near-infrared spe ctra with the TripleSpec 
instrument (jHerter et al.ll2QQ^ on the Apache Point Ob¬ 
servatory 2.5 m telescope to supplement the material de¬ 
scribed above and to aid in the determination of the stel- 


Table 4 

Spectroscopic observations (Keck/HIRES). 


Star 

BJD 

(2,400,000+) 

Exposure 

(sec) 

SNR 

(pix-i) 

RV 

(km s“^) 

KOI-0571 

56135.II242 

900 

18 

-61.7 

KOI-1422 

56192.77304 

900 

13 

-24.8 

KOI-2529 

56472.07615 

900 

13 

-12.2 

KOI-3255 

56472.06602 

729 

24 

-15.3 

KOI-3284 

56450.88292 

2700 

45 

-90.9 

KOI-4005 

56472.05630 

785 

25 

-41.0 

KOI-4087 

56472.04532 

900 

15 

-20.1 

KOI-4427 

56472.03348 

900 

12 

-II.6 

KOI-4622 

56472.08746 

900 

15 

-98.6 

KOI-4742 

56472.02224 

900 

20 

-76.2 

KOI-4745 

56472.09849 

900 

14 

-24.0 


lar properties (Sect. [5j). For both observations, taken on 
UT 2014 June 26, we used the l7l slit giving a resolv¬ 
ing power of R ^ 3500. We alternated exposures in an 
‘ABBA’ fashion, integrating for 4 minutes at each po¬ 
sition until we achieved a SNR of at least 80 per pixel 
throughout the H band (our SNR was slightly higher, 
about 100, for KOI-4622). Our total integrations were 
48 minutes and 64 minutes for KOI-4087 and KOI-4622, 
respectively. We took care to gather spectra of the AO V 
telluric standard star 51Dra (using an exposure time of 
10 s per image in the same ABBA nod strategy) at in¬ 
tervals separated by no more than 30 minutes in time 
from the KOI observations, corresponding to changes in 
airmass of 0.1 from the positions of the KOIs. We also 
gathered flat field images using the bright quartz lamp at 
the Apac he Point Observator y. We used the spextool 
pipeline (jCushing et al.l l2004l ) to co-add flat fields and 
science frames and to extract the spectra from the full- 
frame images. The xt ell cor package was then employed 
to compare the spectra of our standard AO V star 51 Dra 
to that of Vega, in order to identify and remove telluric 
absorption lines in both KOI spectra. 

5. STELLAR PROPERTIES 

We used our Keck/HIRES spectra to estimate the 
effective temperature (Teff), surface gravity (log^), 
and metallicity ([Fe/H]) of our stars using a proce¬ 
dure referred to as SpecMatch (Petigura et ah, in 
prep.). SpecMatch compares an observed target spec¬ 
trum against a set of approximately 800 library spec¬ 
tra obtained with the same instrument that span a wide 
range of parameters (Teff = 3500-7500K, log^ = 2.0- 
5.0 (cgs), and [Fe/H] values from —1.0 to +0.5dex). For 
each comparison SpecMatch places the target spectrum 
onto the same wavelength scale as the library spectrum 
and computes , the sum of the squares of the pixel-by¬ 
pixel differences in normalized intensity. We adopted the 
mean Teff, log^, and [Fe/H] of the ten closest matching 
library spectra (weighted by the y^) as the final set of 
stellar parameters for each target spectrum. They are 
reported in Table [5] for about half of our sample. We 
estimate these parameters to be accurate to about 100 K 
in Teff, 0.10 dex in log^, and 0.10 dex in [Fe/H] based 
on comparisons with benchmark stars having parameters 
measured using detailed LTE modeling. 

For stars cooler than about 4300 K this technique be¬ 
comes less reliable. In that temperature range the library 
of reference spectra is quite sparse, and it is possible for 
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Table 5 

Spectroscopic parameters. 


Star 

^eff 

(K) 

logg 

(cms-2) 

[Fe/H] 

(dex) 

Source 

KOI-0571 

3755 ± 90 


-0.26 ±0.12 

1,2 

KOI-1422 

3572 ± 80 


-0.12 ±0.12 

1,2 

KOI-2529 

4651 ± 100 

4.64 ±0.10 

±0.01 ±0.10 

3 

KOI-3255 

4551 ± 100 

4.67 ±0.10 

0.00 ±0.10 

3 

KOI-3284 

3748 ± 112 


±0.16 ±0.14 

2 

KOI-4005 

5431 ± 100 

4.50 ±0.10 

±0.02 ±0.10 

3 

KOI-4087 

4134 ± 154 


-0.30 ±0.15 

4 

KOI-4427 

3813 ±112 


-0.07 ±0.14 

2 

KOI-4622 

4340 ± 177 


-0.57 ±0.18 

4 

KOI-4742 

4402 ± 100 

4.71 ±0.10 

-0.37 ±0.10 

3 

KOI-4745 

4723 ± 100 

4.62 ±0.10 

-0.01 ±0.10 

3 


Note. — Source s are: 1. IMann et al.l (|2013ai rbh: 2. 

IMuirhead et al.l (|2014l ): 3. SpecMatch; 4. New estimates from 
this paper. 


the SpecMatch Teff value to be biased high simply be¬ 
cause the best matches are hotter. Similar biases are 
likely in log g and [Fe/H]13 Six stars in our sample are in 
this regime: KOI-0571, 1422, 3284, 4087, 4427, and 4622. 
Spectroscopic parameters (except for log g) for the first 
three a nd for K OI-4427 were reported by iMann et al.l 
(|2Q13al lH) and/or iMuirhead et al.f (j2Ql4lFl based on the 
measurement of temperature- and metallicity-sensitive 
features in near infrared i7-band or K-hand spectra. We 
adopted averages of those determinations here. 

For KOI-4087 and KOI-4622 we derived new param¬ 
eters based on the NIR spectra described in Sect. 14.41 
Effective temperatures we re inferred fro n i the empiri¬ 
cal i7-band calibrations of iNewton et al.l (|2Q15f ). which 
were established using stars with bolometric luminosi- 
ties and temperatures from long-baseline interferometry 
(jBovahan et al.l 1201 3[ ) . They are valid between 3100 K 
and 4800 K, and have a typical scatter of 72 K. Spec¬ 
tral indices based on the Mg I doublet at 1.50/im 
and the Al I doublet at 1 .67 gm were measured as de¬ 
scribed by INewton et al.l (j2015[ ). and temperature un¬ 
certainties were obtained from numerical simulations 
following the same work. Metallicities for KOI-4087 
and KOI-4 622 were de t ermined using the NIR cali¬ 
brations of IMann et~ld ] (jM13), which have the ad- 
vanta ge over other calibrati c )ns developed in recent 
years (Roias-Avala et al.l[2010L 120121: iTerrien et al.ll20ll : 
INewton et al.l[2014al) that they were established includ¬ 
ing K dwarfs with effective temperatures greater than 
4000 K. They are thus the most appropriate calibrations 
for t hese two KO I s, whic h lie on the K/M-dwarf bound¬ 
ary. IMann et al.l (|2013c[ ) showed that the TC-band rela¬ 
tions display the strongest correlation with [Fe/H]. We 
adopted their calibration in Eq.(16), which has a scatter 
of 0.11 dex. The analysis of the spectra, the measure¬ 
ment of the spectral features (equivalent widths of the 
2.21 fim Na I doublet, the CO (2-0) band-head in the 
K ba nd, and the H 2 O-K 2 index from lR.oias-Avala et al.l 
[ 20 TI I. as well as the uncertainty estimation followed the 

Even though SpecMatch has difficulty obtaining accurate pa¬ 
rameters for cool stars, it is still easily capable of distinguishing 
dwarfs from giants as the pressure-sensitive lines such as the Mg I b 
triplet are very much narrower in giants than in our dwarf library 
spectra. All of the targets in thi s paper are verified to be dwarfs. 

The IMuirhea^j^^y (|201J) values supersede earlier ones re¬ 
ported by IMuirhead et al.l (|2012l f that used the same spectra. 


procedures described recently bv IMuirhead et ahl (j2014f ). 
who applied the same methodology to a similar set of 
spectra of about 100 M-dwarf KOIs. The new tempera¬ 
tures and metallicities for KOI-4087 and KOI-4622 may 
be found in Table p 

The spectroscopic parameters for our targets were used 
to estimate the stellar properties (primarily the mass and 
radius) by appeali ng to stellar evolut ion models from the 
Dartmouth series (jPotter et al.ll200'^ . follo wing a Monte 
Carlo procedure similar to that described bv lTorres et al.l 
(|2008[ ). While this is straightforward for the stars with 
complete information (Teff, [Ee/H], and log^), we lack a 
log g estimate for the six cooler stars because the spectro¬ 
scopic techniques applied to them do not constrain that 
property. Thus, we cannot establish their precise ages, 
or equivalently, their sizes. This uncertainty is relatively 
unimportant, however, because the radii of cool main- 
sequence stars change little with age. Most of those stars 
show periodic brightness variations that may be inter¬ 
preted as rotational modulation, and in principle knowl¬ 
edge of the rotation periods (Prot) enables one to infer 
a rough age using gyrochronology relations. The age, in 
turn, may be used in place of log^f as a constraint for 
the stellar evolution modeling, although in practice the 
constraint is weak because the radius is not very sensi¬ 
tive to age, as just stated. Our determination of Prot for 
KOI-3284, 4087, 4427, and 4622 is explained in the fol¬ 
lowing section; for KOI-0571 and KOI-1422 we adopted 
the periods measured previously by others (see below). 
O ur age estimates relied on the gy rochronology relations 
of lEpstein fc Pinsonneaultl (|2Q14 i. which are claimed to 
provide values good to ^1 Gyr for stars with masses 
above 0.55 Mq, degrading rapidly below that. The B — V 
color indices required by these calibratio ns were taken 
from t he UBV survey of the Kepler field bv lEverett et al.l 
([ 20 TI I. along with reddening values from the KIC. Gen¬ 
erous uncertainties of 0.1 mag were assigned to the de- 
reddened color indices to account for possible systematic 
errors. Eor KOI-1422 and KOI-3284 the colors were ad¬ 
ditionally corrected for the presence of the close compan¬ 
ions described in Sect. mu assuming they are physically 
associated. The formal ages inferred for KOI-0571, 1422, 
3284, 4087, and 4622 are S.bltlH, 4.3611^;^?, 

1.20^q' 52, and 1 . 82 ^g' 3 ^ Gyr, respectively. As described 
below, KOI-4427 does not display a clear signature of 
rotation, so we have conservatively adopted a broad in¬ 
terval of possible rotation periods (10-45 days) based on 
the full range of periods ob served for M star s of si milar 
temperature, as reported bv lMcOuillan et all (|2013[ ) (see 
their Eig. 10). The corre sponding age range from the gy - 
rochronology relations of lEpstein fc Pinsonneaultl (|2014l ) 
is 1-6 Gyr, which we then used as a constraint for the 
isochrone modeling. 

KOI-0571 and KOI-1422 are systems with multiple 
transiting planet candidates (five each). In such cases, if 
the candidates can be assumed to transit the same star 
in near-circular orbits (required for stability), the transit 
light curve modeling can provide a much stronger con¬ 
straint on the mean stellar density (p^) than for single¬ 
planet candidates because the density estimates from the 
different candidates wit hin the same system can be aver¬ 
aged together (see Sect. 18.ip . Eor these two KOIs we have 
therefore made use of the photometric in our Monte 
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Carlo procedure simultaneously with the age constraint 
to strengthen the determination of the stellar character¬ 
istics. 

The final properties from our stellar evolution mod¬ 
eling are listed in Table [6l computed from the mode 
of the corresponding posterior distributions. Uncertain¬ 
ties correspond to the 68.3% (la) credible intervals from 
the same distributions. Distances were computed from 
the absolute Kg magnitudes inferred from the isochrones 
along with the apparent brig htness from the T wo Mi¬ 
cron All-Sky Survey (2MASS; iCntri et~aI]l2QQ3( ) in the 
same passband, which is the one least affected by in¬ 
terstellar extinction. We nevertheless made corrections 
for extinction that we inferred from reddening values for 
each KOI adopted from the KIC, as above. For the four 
targets with close companions (KOI-1422, 2529, 3255, 
and 3284) the apparent magnitudes were corrected for 
the light contribution of the neighbors. On the assump¬ 
tion that those stars are physically bound to the targets, 
their typical orbital semimajor axes would be approxi¬ 
mately 52, 248, 75, and 64 AU, and their orbital periods 
roughly 550, 4600, 770, and 690 yr, respectively. 

5.1. Stellar rotational periods 

Here we give the details of our determination of the ro¬ 
tation periods used above for four of the cooler targets in 
our sample that have no spectroscopic estimate of log^: 
KOI-3284, 4087, 4427, and 4622. For this application 
the processing of the raw photometry was somewhat dif¬ 
ferent than that described earlier, because we wished to 
retain the astrophysical variations present in the light 
curve while at the same time removing the instrumental 
effects. We therefore used the Presearch Data Condi¬ 
tioning Maximum A-Poste riori (PDC-MAP) data from 
Kepler (jSmith et al.l l201^ . which is designed to meet 
those goals. 

Since the data are unevenly sampled and each quar¬ 
ter has a unique offset, we elected to use a Lomb-Scargle 
style periodogramEU The light curve model is a simple 
sinusoid and thus is linear with respect to the model 
parameters for any trial rotation period, Prot- Using 
weighted linear least squares we are guaranteed to find 
the global maximum likelihood solution at each trial Prot- 
We scanned in frequency space from twice the cadence 
up to twice the total baseline of observations, taking 10^ 
uniform steps in frequency. At each realization, we de¬ 
fined the “power” as (P/Cnuii~P7(7triai)/P^Cnuin where 
BIC is the Bayesian Information Criterion and “null” 
and “trial” refer to the two models under comparison. 
We also performed a second periodogram analysis with a 
finer grid step around any prominent peaks found in the 
original periodograms. 

The resulting periodograms for KOI-3284, 4087, 4427 
and 4622 are shown in Figure IH KOI-3284 and KOI- 
4087 show clear uni-modal peaks at 37.53 ± 0.73 d and 

Stars with active regions have a non-uniform surface bright¬ 
ness d istribution, lea ding to brightness variations as the star ro¬ 
tates (jBuddind 119771b These active regions tend to evolve in lo¬ 
cation and amplitude over timescales of days to years, which can 
cause the periodicities to change as well due to differential rotation 
(| Reinhold et al.l [201311 . Despite the complex nature of individual 
spots, the ensemble population tends to imprint the rotation period 
as a dominant peak in the Fourier domain, allowi ng for an estimate 
of the rotation period using photometry alone (IBasri et al.ll201ll : 
INielsen fc K^75a[2071f . 
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Figure 4. Periodograms of the Kepler PDC-MAP photometry 
for four of the KOIs for which a spectroscopic surface gravity con¬ 
straint is lacking. The insets show a close-up of the main peak, with 
an indication of the full width at half maximum that we assign as 
the error in Prot- No clear result is obtained for KOl-4427. 

17.62 ± 0.08 d, respectively, which we attribute to their 
rotation periods. Here we have taken the la confi¬ 
dence intervals to be the full width of the peak in the 
periodogr am at half maximum. For comparison, we 
note that iMcQuillan et al.l (j20I3f ) used an autocorrela¬ 
tion periodogram and found peaks at similar periods 
(37.026d and 17.54 2 d). Addition a lly, K OI-4087 is in the 
sample studied by INielsen et al.l (j20I3f ). who obtained 
Prot = 17.496 d, in good agreement with the other re¬ 
sults. KOI-4427 shows no clear peaks above the noise 
l evel, in agreement with the findings of IMcQuillan et ^ 
so we were unable to assign a rotation period in 
this case. A possible reason for complex structure of the 
periodogram in this case is the presence of other unre¬ 
solved stars in the aperture. 

KOI-462 2 exhibits two p e aks a nd is no t included in 
either the IMcQuillan et al.l (|2QI3D or the INielsen et al.l 
(|2QI3D samples. The dominant peak occurs at 21.85 ± 
0.22 d, with a second peak at ^25% lower power with 
P = 10.72 ± 0.06 d. We note that this second peak oc¬ 
curs not only at a lower power but also at nearly exactly 
one half the period of the higher peak. We argue that 
the longer period is that of the stellar rotation, and the 
shorter period is an alias caused by more than one active 
region on the stellar surface. 

Two other cool KQIs have had their rotation periods 
determined previously (33.7 0 T Q.3I d for KOI-Q57 I and 
36.5 ± 17.6 d for KQI-1422: IMcQuillan et aTll2QI3D . and 
we adopt those values here as published. 


6. CANDIDATE VALIDATION 

The Doppler signals (velocity semi-amplitudes) ex¬ 
pected of the candidates in our sample, if due to planets 
around the targets, range from 0.4 to 1.4 ms“^ based 
on rough estimates of their masses using the prelim- 
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Table 6 

Stellar properties. 



Age 



logs 

pi. 


Mv 

Mk, 

Distance 

Star 

(Gyr) 

(Me) 

(Rq) 

(cm s“^) 

(gcm-3) 

+ 0 ) 

(mag) 

(mag) 

(pc) 

KOI-0571 

4 Q+O-6 

0 544+^-^2'^ 
U.O^^_Q Q21 

n c;9Q+0.023 

4.736l°;°20 

"^•^^-0.39 

0 055+0-011 

U.UOO_Q QQg 

^•01-0.24 

5 41+0.14 

O-^l-0.15 

172+13 

I'^-lO 

KOI-1422 

^•^-1.6 

0 454+0-033 

U.^ZD_o 027 

4.8331^25 

7 94+1-3'^ 

0 027+®-®®® 

U.UZ/ _Q QQ4 


6 07+0-20 
O-O' -0.23 

226t.1l 

KOI-2529 

3 

3 

u. q29 

0.697l°;°28 

4:.Diy_Q q28 

3 03+0-24 
o.uo_0 25 

0 199+®-®^^ 
u.iyy_Q q25 

6.98l°;25 


618134 

KOI-3255 

2 9+'^-^ 

^•^-0.3 

0 707+0-033 

O'- ' "O' _0.027 

0.6801°;°26 

^.u^a_Q Q2g 

q 1 q+0.25 
'^•-*-^-0.25 

n 1 7Q+0.035 

U.i /O_o 022 

/ .ZU_o 24 


4171 

KOI-3284 

4 4+0-8 


0.52o1°;“8 

4 740+0-059 
<^u_o 029 

5 59+1.53 

0-0^_0.77 

0 044+0-017 

O-O44_0 012 

9 55+0.54 
y.00-0 44 

O.OU_o 25 

1451^^ 

KOI-4005 

7 0-1-3.6 
‘•^-3.9 

0-884l°;“| 

0.866l°;°I8 

4 51 4+0-035 
4.0i4_Q Q73 

i.o»_o 38 

0 5ei +0.153 
U.DOi_o 079 

^•^0_o.25 


6931“ 

KOI-4087 

1 q-|-0.6 

0 575+0.043 
U.O/0_Q Q47 

0 550+0.029 
U.009_q Q 54 

4 706+0-049 

' ’0O_0.016 

4 76+^-®^ 
'0-0.48 

0 079+®-®^^ 

a_o.022 

0 qq + 0.55 
O.OO-o 32 

o.uz_o 13 

26111^ 

KOI-4427 

3'6lt3 

'4.0ZD_q Qg2 

0 505+0.038 
'4.0UO_q Qgg 

4 751+0-067 

/OJ._o 030 

5 79+1-87 

o- ' ^-0.82 

0.043ini2 

9.34l“l[ 

5 54+0-42 

0-0"1-0.24 

240132 

KOI-4622 

i-stoi 

n 579+0.049 
U.O/Z_o 053 

0.5501°;“8 

4 7. .+0.047 

4. /'i0_o 024 


0.089l°;°“ 

8-02l!^:^l 

O.UU-0 20 

284128 

KOI-4742 

^.y_0 2 

'4.DUy_Q q26 

0.598l°;“3 

"±-O<o_0.021 

4 01+0-37 
^•01-0.30 

0 117+0-024 
0-11'-0.016 

7 7q+0.28 

' •'O_o.25 

4 79+°-^^ 

342l« 

KOI-4745 

q 9+7.5 

n 7qo+0.033 

U. /00_Q Q2g 

0.706l°;“8 

^.UJ.^_0 029 

^.au_Q 25 

n 91 7+0.043 
‘ -0.027 

6.83l°;i 

^•'^^-0.11 

779l« 


Note. — Ages, log^f, and p-k are best-fit values from the Dartmouth models, constrained by either a spectroscopic log^f, a 
light-curve-derived and/or a gyrochronology age based on the rotation period (see text). 


inary radius measures and the mass-radius relation of 
IWeiss fc Marcvl (j2Q14l ). Given the faintness of the par¬ 
ent stars {Kp = 14.3-15.9), the detection of such small 
variations would be very challenging, and therefore the 
planetary nature of these objects cannot presently be 
“confirmed” in the usual way, i.e., by establishing that 
the orbiting objects are of planetary mass. Instead they 
must be “validated” statistically, by showing that the 
likelihood of a true planet (which we refer to hereafter 
as the ‘planet prior’) is orders of magnitude larger than 
that of a false positive. We describe this process below. 

6.1. Overview 

The types of astrophysical false positives we consider 
in our analysis involve other unseen stars in the photo¬ 
metric aperture of Kepler that are eclipsed by an orbiting 
object, and have their eclipses attenuated by the light of 
the target such that they mimic shallow planetary tran¬ 
sits. Examples of these situations include background 
or foreground eclipsing binaries (‘BEB’), background or 
foreground stars transited by a (larger) planet (‘BP’), 
and physically associated stars transited by a smaller 
star or by a planet. Such physically associated stars 
will typically be close enough to the target that they 
are generally unresolved in high-resolution imaging. We 
refer to these hierarchical triple configurations as ‘HTS’ 
or ‘HTP’, depending on whether the object orbiting the 
physical companion is a star or a planet 0 

The procedure we use d for validation is BLENDER 
([Torres et all I2QQ4 I2Q11I: lEressin et al.l |2Q12[) . which 

We point out that other statistical studies using Kepler data 
(e.g., IMorton &: JohnsQnll20Hl : IMQrtQnll2012l : ILissauer et al.|[2014l) 

have not considered the HTP scenario as a false positive, arguing 
that the exact location of the planet (whether on the intended tar¬ 
get or a bound companion) is inconsequential although its size may 
be larger if around the companion. Validation without consider¬ 
ing these kinds of false positives is significantly easier, particularly 
since they tend to dominate the blend frequency, as we show be¬ 
low. For the present work, however, planet size is critical because 
it affects habitability: a planet that transits an unseen companion 
instead of the target may be too large to be rocky, and is thus les s 
interesting for our purposes (see, e.g.. IAlibertll2014 : IRogersll2014l) . 
We therefore count HTP configurations as blends. 


has been applied successfully in the past to many 
of the most i nteresting candidates revealed by Ke¬ 
pler (see, e.g., [Borucki et al.l 12013 : iBarclav et al.l l2Q13l: 
Meibom et al.l I2Q13I: iBallard et al.l I2Q13I: iKipping et al l 


2Q14b[ ). BLENDER makes full use of the detailed shape 


of the transits to limit the pool of viable blends. It does 
this by simulating large numbers of blend scenarios and 
comparing each of them with the Kepler photometry in 
a sense. Eits that give the wrong shape for the tran¬ 
sit are considered to be ruled out. This enables us to 
place useful constraints on the properties of the objects 
that make up the blend, including their sizes and masses, 
overall color and brightness, the linear distance between 
the background/foreground eclipsing pair and the tar¬ 
get, and even the eccentricities of the orbits. Those 
constraints are then used to estimate the frequencies of 
blends of different kinds. We note that the simulated 
light curves generated by BLENDER take full account of 
any known extra light in the aperture, such as that com¬ 
ing from the close companions reported earlier (see Ta¬ 
ble [2]). We also point out that here we only consider 
blends involving main-sequence stars, as BEBs with a 
giant component produce light curves with a shape that 
does not mimic a true planetary transit when observed 
with such high precision as delivered by Kepler. 

The photometric data we use here are the long-cadence 
PDC-MAP time series, as in Sect. EH detrended to 
remove signals at time scales long compared to the 
transit durations. Throughout this analysis we fol¬ 
lowed the nomenclature established in previous BLENDER 
studies, designating the objects in the eclipsing pair 
as the “secondary” and “tertiary”, and the target it¬ 
self as the “primary”. We drew stellar properties for 
the primary, secondary, and tertiary (masses, radii, 
and absolute brightness in the Kepler and other pass- 
band s) from model isoc hrones from the Dartmouth se¬ 
ries ([Potter et al.l [200^ : the selection of the isochrone 
for the primary was based on the spectroscopic proper¬ 
ties given earlier. The BLENDER studies cited above may 
be consulted for further technical details. 

In summary, our validations proceed in two stages. 
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First, we use BLENDER proper to derive constraints on 
blend scenarios from their light curve shapes. Then, 
we incorporate other constraints from follow-up obser¬ 
vations to compute the blend frequencies and the planet 
prior by means of Monte Carlo simulations. We now 
address each of these steps in turn. 

6.2. BLENDER constraints 

To illustrate the constraints derived from the shape of 
the transit, we focus here on one of our candidates, KOI- 
4005.01. BLENDER indicates that background eclipsing 
binaries are only able to produce viable false positives if 
the primary star of the eclipsing pair is restricted to a 
narrow range of main-sequence masseS between about 
0.7 Mq and 1.4 Mq, and to a limited interval in bright¬ 
ness {Kp magnitude) relative to the target correspond¬ 
ing to AKp < 5.5. Figure [5] shows the landscape 
for all blends of this kind (BEB scenario) in a repre¬ 
sentative cross-section of parameter space. Regions out¬ 
side of the 3(7 contour correspond to configurations with 
light curves giving a poor fit to the observations, i.e., 
much worse than a true planet fit. These blends are 
therefore excluded. The figure also illustrates some of 
the additional constraints available from our follow-up 
observations for this candidate. For example, analysis 
of our Keck/HIRES spectra of KOI-4005 generally rules 
out companions within 5 magnitudes of the primary if 
they are angularly close enough to fall within the spec¬ 
trograph slit 13 This can eliminate much o f p arameter 
space (see green hatched area in the figur43)- Addi¬ 
tionally, by comparing the r — Kg colors of the simu¬ 
lated blends with the m easured color index of the target 
{r-K, = 1.58n + n.n29: IRrnwn et, alllaTllI) . we find that 
some of the BEB scenarios we have simulated are either 
too blue or too red by more than 3(7 (blue hatched areas 
in the figure), and are therefore also excluded. How¬ 
ever, for this particular KOI all BEB blends with the 
wrong color are already excluded by BLENDER for giv- 
ing poor fits (i . e., the y are outside of the 3(7 contour). 
iSanterne et al.l (|2Q13[ ) have pointed out that an addi¬ 
tional source of blends involves eclipsing binaries with 
eccentric orbits that are oriented so that they show only 
a secondary eclipse, as viewed by Kepler. More gener¬ 
ally, we note that they could also show only a primary 
eclipse. We find, though, that while the depth and shape 
of the diluted eclipses may indeed match the transit sig¬ 
nal in some cases, the combined color and brightness of 
these blends are such that they are generally ruled out by 
the spectroscopic constraint and/or the measured color 
of the KOI. 

Stars more massive than about 1.4 Mq will generally have 
evolved to become giants for ages typical of the field (~3 Gyr), and 
such stars yield light curves that have a very different shape than 
a transit, as stated earlier, so they do not constitute viable blends. 

For typical seeing conditions at the Keck telescope stars that 
are beyond 0743 can still imprint their lines on the target spec¬ 
trum, and be detected, as can stars at wider separations that hap¬ 
pen to be aligned along the slit. Nevertheless, to be conservative 
we assume here that these companions are spectroscopically unde¬ 
tectable. 

Note, however, that not all blends in the green hatched area 
are excluded. This is only the case if the background star falls on 
the slit, and its RV is more than 10kms“^ different from that of 
the target. In all other cases we consider the blends to be viable, 
even if they have AKp < 5. 
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Figure 5. Map of the surface (goodness of fit) for KOI- 
4005.01 corresponding to blends involving background eclipsing 
binaries. On the vertical axis we represent the linear distance 
between the BEB and the target (Dbeb — T^targ), expressed 
for convenience in terms of the difference in distance modulus, 
A(5 = 5 log(DBEB/L)targ)- Only blends within the solid white con¬ 
tour (darker colors) provide fits to the Kepler light curve that are 
within acceptable limits (3cr, where a is the significance level of 
the x^ difference compared to a transiting planet model fit; see 
IFressin et al.ll201^) . Other concentric colored areas (lighter colors) 
represent fits that are increasingly worse (4cr, 5cr, etc.), which we 
consider to be ruled out. The blue cross-hatched areas correspond 
to regions of parameter space where the blends are either too red 
(left) or too blue (right) compared to the measured r — Kg color 
of the target, by more than three times the measurement uncer¬ 
tainty. The dashed green line labeled AKp = 5.5 is tangent to the 
white contour from above and corresponds to the faintest viable 
blends. The green line labeled AKp = 5.0 represents the spectro¬ 
scopic limit on faint background stars. All simulated blends below 
this line (green hatched region) are brighter and are generally ex¬ 
cluded if the BEB is angularly close enough to the target to fall 
within the slit of the spectrograph. Thus, very few blends remain 
viable. 


Eor blends involving a background or foreground star 
transited by a planet (BP scenario) there is a wide range 
of secondary masses that produce good fits to the Kepler 
photometry of KOI-4005, as shown in Eigure[6l BLENDER 
indicates that the faintest of these blends are about 5.8 
magnitudes fainter than the target in the Kp passband. 
In this case, however, the spectroscopic and color con¬ 
straints drastically reduce the pool of viable false posi¬ 
tives. 

The x^ map for blends involving a physically associ¬ 
ated companion to KOI-4005 transited by a larger planet 
(HTP scenario) is seen in Eigure [71 and shows the size 
of the tertiary as a function of the mass of the com¬ 
panion star (secondary). In this case BLENDER restricts 
the false positives to be in a narrow strip of parame¬ 
ter space corresponding to secondary masses larger than 
about 0.25 Mq, and planetary sizes between 0.25 and 
about O.SRjup (2.8-9.0 Rq). As in the other scenarios, 
color and brightness constraints allow us to reject many 
of these blends. 

Einally, BLENDER indicates that pairs of eclipsing stars 
orbiting the target (HTS scenario) invariably produce 
light curves with the wrong shape for a transit, or fea¬ 
ture noticeable secondary eclipses that are not observed 
in the photometry of KOI-4005.01, or, if they show only 
a single ecl ipse clue to a high e ccentricity and special 
orientation (jSanterne et al.|[2QT^ . the overall brightness 
would make the eclipsing binary detectable and/or its 
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Figure 6. Similar to Figure [S] (and with the same color scheme) 
for blends involving background or foreground stars transited by 
a planet (BP scenario). The faintest blends giving acceptable fits 
have AKp = 5.8 relative to the target (dashed green line). Blends 
below the solid green (AKp = 5.0, hatched green area) are generally 
excluded by the spectroscopic constraint unless the intruding star 
is more than 0''43 from the target, or closer than 0^43 but with 
ARV < 10kms“^ (see text). 
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Figure 7. Similar to Figure [5] for the case of physical companions 
to KOI-4005 that are transited by a planet (HTP). Only companion 
stars with masses larger than about 0.85 Mq or smaller than about 
0.42 Mq yield blend colors that are consistent with the measured 
r — Kg index of KOI-4005. However, many of these blends are 
eliminated by the spectroscopic constraint, as they are brighter 
than AKp = 5.0. 

color inconsistent with the measurement. These types 
of configurations are therefore easily excluded. A simi¬ 
lar result was found for each of the other candidates in 
our sample, and indeed all previous BLENDER analyses of 
KOIs have also found that HTS blends are always ruled 
out with photometry of the quality delivered by the Ke¬ 
pler instrument, when combined with the observational 
constraints. Consequently, we do not consider HTS sce¬ 
narios any further. 

Maps of the landscapes for the BEB, BP, and HTP 
scenarios analogous to those described above, along with 
the additional restrictions based on the color and the 
spectroscopic brightness limits, are shown for each of our 
other candidates in Figure set ISHTOl 

6.3. Estimating blend frequeneies and the planet prior 


The BLENDER constraints just described significantly 
reduce the number of viable BEB, BP, and HTP blends, 
but not all of the ones that remain can be ruled out by 
our follow-up observations. The next step of the analysis 
is therefore to quantify their expected frequencies, which 
we did via Monte Carlo simulations. Specifically, we gen¬ 
erated large numbers of realistic false positives according 
to each of the blend scenarios, we eliminated those that 
BLENDER tells us would not match the Kepler light curve 
or that would have been detected in our follow-up obser¬ 
vations, and we counted the survivors. These numerical 
experiments rely on the known distributions of binary 
star properties, the number density of stars near each 
target, and estimates of the rates of occurrence of tran¬ 
siting planets and of eclipsing binaries from the Kepler 
Mission itself. Those rates of occurrence (as well as any 
dependence they may have on orbital period or other 
properties) are implicit in the lists of KOI s and eclips¬ 
ing b i naries generated by the Kepler team (jBurke et al.l 
I2Q14I: iSlawson et al.ll2QTlh that we used in our analysis, 
when normalized by the total number of targets observed 
by the spacecraft. Since we have recently made sig¬ 
nificant enhancements in these Monte Carlo procedures 
compared to previous applications, we take the opportu¬ 
nity to describe them in some detail here. 

For the HTP scenario we simulated companion stars to 
each target following the distributions of binary proper¬ 
ties (mnssj;atioSj_or^bita][^riods , eccentrichies)_m^ 
by iRaghavan et al.l (j2QlQ[ ) and iTokovininI (120141 ). We 
placed the companions in random orbits around the tar¬ 
get assuming isotropically distributed inclination angles 
(i) and uniformly distributed longitudes of periastron 
(cj), as well as random orbital phases. Other relevant 
properties of these simulated companions (size, absolute 
magnitudes, and colors corresponding to their simulated 
masses) were inferred from the same isochrone used for 
the target, as the two stars are assumed to be coeval in 
this scenario. In particular, we computed the distance 
to the system by imposing the condition that the to¬ 
tal apparent brightness of the two stars (accounting for 
extinction) must be the same as the Kp magnitude as¬ 
sociated with the KOI, and this in turn enables us to 
compute the angular separation. We then assigned to 
each of these secondary stars a random transiting planet 
drawn from the actual list of KOIs hosted at the NASA 
Exoplanet Archive (list downloaded on 2014 March 26), 
which ensures that their properties (including any cor¬ 
relations among them) are as realistic as possible0 We 
accepted only planets with periods similar to the KOI 
under investigation (within a factor of two). The ratio¬ 
nale for this is that the relevant HTP blend frequency 
is that of configurations that involve planets with peri¬ 
ods near that of the candidate, as those frequencies de¬ 
pend strongly on orbital period. We then examined the 
properties of the stellar companions and their planets, 
and rejected scenarios that do not satisfy the BLENDER 
restrictions on companion mass, planetary size, and or¬ 
bital eccentricity of the planet (related to the transit 

As pointed out by ILissauer et al ] lEniD, when attempting 
to validate targets with a single transiting planet candidate it is 
more appropriate to exclude from the list of known KOIs all stars 
with multiple candidates, as this unduly inflates the rates of occur¬ 
rence of planets (particularly small ones that are more common in 
multiple systems). We have followed this recommendation here. 




























Validation of small HZ planets 


13 


0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Mass of background star (Mq) 






0.2 0.4 0.6 0.8 1.0 

Mass of companion star (Mq) 


0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Moss of background star (Mq) 


Moss of background star (Mq) 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Moss of background star (Mq) 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Moss of background star (Mq) 


0.2 0.4 0.6 0.8 1.0 

Mass of companion star (Mq) 



0.4 0.6 0.8 1.0 

Mass of background star (Mq) 


0.2 0.4 0.6 0.8 1.0 

Mass of companion star (Mq) 


Background eclipsing binaries 

- KOI-2529.02 


0 

0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Mass of background star (Mq) 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Mass of background star (Mq) 


0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Mass of background star (Mq) 


0.2 0.4 0.6 0.8 1.0 

Mass of companion star (Mq) 


Figure 8. BLENDER landscapes for five KOI targets. Each row corresponds to the KOI labeled in the left panel, and shows the three 
blend scenarios. Observational constraints from spectroscopy and color are also shown. See Figures [SHZl for a detailed description. 
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Figure 9. Same as Figure [8]for a second group of five KOIs in the sample. 
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Figure 10. Same as Figure [8] for the remaining KOI in the sample. 

duration), as they yield poor fits to the transit light 
curve. We also rejected configurations with companions 
bright enough that they would have been detected from 
our high-resolution imaging observations, our spectro¬ 
scopic observations, by our centroid motion analysis, or 
that would yield an overall color for the bl end inconsis¬ 
tent with the measured color of the target (|Brown et al.l 
unni). In applying the spectroscopic constraint we ex¬ 
cluded companions brighter than the limits described in 
Sect. 14.31 only if they are within 0'.'43 of the target (half 
width of the spectrometer slit), unless their radial veloc¬ 
ity computed from the simulated orbit around the target 
is within 10kms“^ of that of the primary. In that case, 
blending of the spectral lines of the secondary with those 
of the primary would prevent detection, so we considered 
those blends still viable. Finally, we retained only the 
false positive configurations that wo uld be dynamically 
stable according to the criterion of iHolman V WiegertI 
(HHi). Since not every stellar companion will necessar¬ 
ily have a transiting planet and therefore act as a blend, 
as we assumed above, we adjusted the blend frequencies 
by accounting for the rate of occurrence of planets of each 
size and period as inferred from the KOI list itself. Minor 
corrections were applied to those rates for incompleteness 
and for the incidence of false positives in the KOI list, 
by me ans of simulations as described by iFressin et al l 
(j2QI3[ ). An additional adjustment was made to account 
for the dependence of the transit probability on the size 
of the secondary relative to the average size of the host 
stars in the KOI list. The final blend frequencies were 
obtained by multiplying by the occurrence rate of binary 
stars wit h planets as a fmmtion of binary semimajor axis, 
following IWang et~al] (j2QI4l ). 

The calculation of the frequency of background or fore¬ 
ground stars transited by a planet (BP) proceeded in 
a similar fashion, and scales with the number density 
of stars in the vicinity of each target, for which we re- 
lied on pred i ctions from the Galactic structure model of 
iRobin et~al] (j2QQ3L l2QI2[ ) (Besancon model) 0 We gen- 

We note that the current online version of the Besangon simu¬ 
lator (available at ht tp: //model. obs- besancon. f r ) relies on stellar 
evolution models bv iHavwoodl (|1994l ) to generate some of the stel¬ 
lar properties returned, including absolute magnitudes, colors, tem¬ 
peratures, and radii. Many of the simulated main-sequence stars of 
interest for our work are significantly coole r than the Sun (spectral 
type K and M), and the IHavwoodl (| 199411 models are not specifi¬ 
cally designed for such stars as they adopt gray boundary condi¬ 
tions and a rather simple equation of state, both of which result in 
incre asingly biased predictions for cool stars fe.g.. [Chabrier et al.l 
[Mnh . Because of this, we have preferred to recompute the relevant 
properties of all simulated stars adopting from the output of the 


crated a list of simulated stars in a 5 square-degree area 
around each of our targets, including their kinematic 
properties (radial velocity), and drew stars randomly 
from this list assigning them a random angular separa¬ 
tion from the target within the corresponding 3cr exclu¬ 
sion region derived from our centroid motion analysis 
(since stars outside this area would have been detected; 
see Sect. 14.2[) . We assigned a random planet to each of 
these stars that we drew from the KOI list, keeping only 
those with orbital periods within a factor of two of that of 
the candidate, as done previously. False positive scenar¬ 
ios that do not meet the constraints from BLENDER were 
rejected, along with those that would have been flagged 
by our high-resolution imaging observations or color in¬ 
formation. Further blends were excluded by the spectro¬ 
scopic limits on the brightness of any unseen companions, 
as before, except for background/foreground stars whose 
spectral lines would be blended with those of the pri¬ 
mary (i.e., those with simulated radial velocities within 
I0kms“^ of the measured heliocentric velocity of the 
target), which we regarded as viable. Surviving blends 
were weighted by the corresponding rate of occurrence of 
transiting planets of that size and period, corrected as in 
the HTP case for incompleteness and false positives and 
adjusted for the dependence of the transit probability on 
the secondary radius. The final blend frequencies were 
then normalized by the ratio of the areas between the 
centroid exclusion region and 5 square degrees. 

To estimate the frequency of BEBs acting as blends 
we again drew stars randomly from a list generated 
with the Besancon Galactic structure model, and paired 
them with an orbiting stellar companion (tertiary) gen- 
erated from the distr ibutio ns of bina r y pro perties by 
iRaghavan et al.l (j2QIQf ) and iTokovininI (j2QI4f ). We as¬ 
signed orbital periods to these secondary-tertiary pairs 
drawn ran domly from th e cata log of Kepler eclipsing 
binaries of iSlawson et al.l (j2QIl[ ) , retaining only those 
within a factor of two of the period of the candidate. 
Each blend was also assigned a random angular separa¬ 
tion from the target within the 3cr exclusion limit from 
our centroid motion analysis. After rejecting scenarios 
that do not meet the constraints from BLENDER for ac¬ 
ceptable fits to the transit light curve, we applied the 

Besangon simulator only the basic characteristics that derive from 
the initial mass function and the Galactic structure model, which 
are the stellar mass, distance, age, and met allicity. We then u sed 
model isochrones from the Dartmouth series (jPotter et al.ll2008l i to 
generate all other stellar quantities for our purposes, as the phys¬ 
ical ingredients of these isochrones give more realistic predictions 
for the lower main-sequence. 
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observational constraints (imaging, spectroscopy, color 
information) as in the BP scenario, and tallied the re¬ 
maining viable blends. For computing the final BEB fre¬ 
quency, each surviving blend was assigned a weight given 
by the ratio of the eclipse probability of the particular 
blend (a function of the secondary and tertiary sizes) 
to the average eclipse probability for the iSlawson et al l 
(ImTh sample as a whole. 

The probability that the candidate is a true planet as 
opposed to a false positive (‘planet prior’, PL) was es¬ 
timated by simply counting the number of known KOIs 
with periods within a factor of two of each candidate, 
radii within 3cr of the measured planetary radius, and 
with similar transit durations. Incompleteness and false 
positive corrections were applied in the same way as 
above. 

6.4. Results 

Blend frequencies for the HTP, BP, and BEB scenar¬ 
ios are presented for each of our candidates in Table [3 
along with the planet priors, PL. Eor the present work 
we have set a threshold for validation equivalent to a 3cr 
confidence level, consistent with previous applications of 
BLENDER. The statistical significance may be expressed 
as PL/(PL + HTP -1- BP + BEB). Eor a 3cr valida¬ 
tion we therefore require a planet prior that is at least 
1/(1/99.73% — 1) « 370 times larger than the total blend 
frequency. We list these ‘odds ratios’ in Table El with 
the corresponding significance levels. Eor eleven of the 
twelve candidates we have achieved very robust valida¬ 
tions at the 3(7 level or higher, strongly supporting their 
true planetary nature. On this basis, the eight that have 
not been previously validated by others are given the new 
planet designations Kepler-436 b to Kepler-443 b (see Ta¬ 
ble [1]). KOI-4427.01 has a lower odds ratio of only ^100, 
corresponding to a significance level slightly over 99%. 
To avoid confusion in the rest of the paper we will con¬ 
tinue to refer to the newly validated planets by their KOI 
numbers. 

Strictly speaking, the above results compare the likeli¬ 
hood that the transit signals originate on the target itself 
with the likelihood that they come from eclipses asso¬ 
ciated with some other unseen star in the photometric 
aperture. However, four of our target stars (KOI-1422, 
2529, 3255, and 3284) have known close companions from 
our high-resolution imaging observations, none of which 
can be ruled out as the source of the t ransit signal by our 
centroid motion analysis in Sect. 14.21 Thus, the possibil¬ 
ity remains that the planets causing the signals orbit the 
companions rather than the targets. In that case, the 
inferred radius of the planets would be somewhat larger 
than the values derived in the next section, the exact 
amount depending on the brightness difference (dilution 
effect) and the physical size of the companion stars. 

An estimate of the likelihood that the planets tran¬ 
sit the companions was derived as follows. In each of 
the four cases we made the assumption that the com¬ 
panions are physically associated with the target. This 
is reasonable, as the probability of having an unrelated 
background star as close to the primary stars as observed 
is very s mall compared to that of a physical companion 
(see also lHorch et al .1120141 ). based on the number density 
of stars around each KOI from the Bes ancon model and 
the known rate of occurrence of binaries (jRaghavan et al.l 


[20T(1 iWang ^nil[2nTl . Eurthermore, no signs of any of 
these companions are seen in our Keck/HIRES spectra 
even though in most cases they are bright enough to have 
been detected. This implies that their radial velocities 
must be similar to those of the primary stars so that the 
spectral lines are blended, preventing detection. This 
again argues for physical associatiouEj In each case we 
then used the same isochrone as for the primaries to infer 
the companion properties, and we computed the ‘planet 
prior’ in the same way as for the targets, but this time 
assuming the planet transits the companion (i.e., increas¬ 
ing the dilution factor by the appropriate amount). We 
represent this planet prior for the companions as FLcomp^ 
and we report these values in Table El Also given in the 
table are the probabilities that the planets transit the 
target rather than the companion, which we computed 
as 7^[targ] = PL/(PL -|- PLcomp)- For KOI-2529.02 the 
results allow us to state with high confidence that the 
planet orbits the target. Eor KOI-I422.04, KOI-I422.05, 
and KOI-3284.0I the calculations indicate a slight or 
modest preference for the planet being around the tar¬ 
get rather than the companion, whereas for KOI-3255.0I 
the result is inconclusive. We return to this issue later 
in Sect. 18.31 

7. LIGHT CURVE FITS AND PLANETARY PARAMETERS 

Having validated eleven of our candidates as bona-fide 
planets with very high confidence, and KOI-4427.01 to 
a lower degree of confidence, we subjected their light 
curves to a detailed modeling to infer the properties 
of the transiting objects. Eor this we used the de¬ 
trended SAP time series from Kepler described earlier 
in Se ct. [3l The tran si ts we re modeled using the stan¬ 
dard iMaiidel fc Agoll (j2QQ2f ) algorithm employing the 
quadratic limb-darkening law. This simple model as¬ 
sumes a spherical, opaque planet transiting a spherically 
symmetric luminous star on a circular Keplerian orbit. 
We re-sampled the long-cadence data onto short -cadence 
sampl ing following the method described by iKippind 
(l2nil . to avoid smearing effects. Our model has eight 
free parameters in total. These are the orbital period, P, 
the time of transit center, r, the planet-to-star radius ra¬ 
tio, Rp/R^, the mean stellar density, p*,phot: the impact 
parameter, 6, the logarithm of fractional light contamina¬ 
tion, log yd, and the quadratic limb-darkening coefficients 
qi and q 2 . All of these parameters have uniform priors in 
our fits, except p*,phot 5 for which we employed a Jeffrey’s 
prior between 10“^ and lO^gcrn”^, and log yd, for which 
we used an informative Gaussian prior where appropri¬ 
ate. 

In our model, yd is the flux of any contaminating sources 
in the aperture divided by the flux of the target star. The 
purpose of this term is to account for dilution of the tran¬ 
sit light curve due to close companions identified by the 
high spatial resolution imaging observations. We com¬ 
puted these ratios here based on the magnitude differ¬ 
ences of the companions listed in Table [21 transformed to 
the Kp band, and assuming that all their flux is included 
in the aperture. Eor the wider companions (4-5" separa¬ 
tions) that are also in the aperture the transformed mag- 

In an independent study of KOI-I422. IStar et al.l (|20I4I I also 
reached the conclusion that the companion is physically bound, 
based on numerical experiments with model isochrones. 
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Table 7 

Blend frequencies, planet priors, odds ratios, and significance level of the validation for our targets. 


Candidate 

HTP 

BP 

BEB 


PL 

Odds ratio 

Significance 

P Lcomp 

P[targ] 

KOI-0571.05 

1.15 

X 

10“® 

1.32 X 10-7 

6.60 X 10-® 

2.15 

X 

10-9 

1595 

99 . 94 % 



KOI-1422.04 

7.75 

X 

10-7 

4.74 X 10-® 

1.28 X 10-® 

3.19 

X 

10-1 

382 

99 . 74 % 

5.94 X 10-1 

78.4% 

KOI-1422.05 

4.49 

X 

10-7 

4.84 X 10-® 

1.29 X 10-® 

4.63 

X 

10-1 

907 

99.89% 

7.25 X 10-® 

86.5% 

KOI-2529.02 

3.24 

X 

10-® 

2.17 X 10-9 

5.55 X 10-9 

1.31 

X 

10-9 

403 

99.75% 

5.61 X 10-® 

99.6% 

KOI-3255.01 

4.29 

X 

10-7 

4.79 X 10-“ 


4.91 

X 

10-1 

1144 

99.91% 

4.39 X 10-1 

52.8% 

KOI-3284.01 

6.71 

X 

10-7 

1.10 X 10-9 

9.27 X 10-19 

1.57 

X 

10-9 

2333 

99.96% 

8.02 X 10-1 

66.2% 

KOI-4005.01 

3.71 

X 

10-7 



6.61 

X 

10-1 

1782 

99.94% 



KOI-4087.01 

1.34 

X 

10-7 



1.32 

X 

10-1 

985 

99.90% 



KOI-4427.01 

2.23 

X 

10-® 

1.07 X 10-9 

1.39 X 10-19 

2.62 

X 

10-1 

117 

99.16% 



KOI-4622.01 

9.98 

X 

10-® 

3.82 X 10-19 

6.68 X 10-12 

2.08 

X 

10-9 

20761 

99.99% 



KOI-4742.01 

1.37 

X 

10-® 

4.89 X 10-19 

5.45 X 10-12 

1.26 

X 

10-9 

919 

99.89% 



KOI-4745.01 

1.53 

X 

10-® 

1.41 X 10-9 

9.19 X 10-12 

2.85 

X 

10-9 

1861 

99.95% 




Note. — For several of the candidates the blend frequencies for the BP and/or BEB scenarios are negligible and are not listed. 
For the four targets with close companions PLcomp is the planet prior assuming the planet transits the close companion, and P[targ] 
the probability that the planet transits the target rather than the companion. 


nitudes were taken from the CFOP Web site, and for the 
closer ones we carried out the conversion using the Dart¬ 
mouth isochrone for each target, on the assumption that 
the companions are physically associated. In the case 
of KOI-1422 we adopte d the magn i tude d ifference in the 
Kp band as reported bv iStar et al.l (j2Q14l ). The P factors 
range from 0.047 ± 0.009 for KOI-2529 to 0.83 dz 0.18 for 
KOI-3255, which has a very bright companion. Addition¬ 
ally, other nearby Kepler sources may contaminate the 
light curve (e.g., bleeding into the photometric aperture) 
and were corrected for using the reported contamination 
factors from the Kepler pipeline, treating them as fixed 
parameters. 

It is worth noting that we did not directly fit for the 
standard quadratic limb darkening coefficients ui and 
1 ^ 2 , but instead we use d the transforme d parameters qi 
and q 2 as advocated bv iKippinj (j2013b[ ). in order to im¬ 
pose efficient, uninformative, and physical priors for the 
limb darkening profile. We note also that the out-of- 
transit baseline flux for each transit epoch was fitted as 
well. However, in this case we used a linear minimiza- 
tio n to derive the baseline flux, similar to that described 
bv iKnndnrthv et al.l (|2Q11[) . This treats the baseline flux 
simply as a nuisance parameter that is not marginalized 
against, but rather minimized at each Monte Carlo re¬ 
alization, reducing the number of free parameters yet at 
the same time allowing us to account for any residual 
offsets in the detrended light curves. 

To regress our 8 -parameter model to the observations, 
we employed the multi-mo dal nested sampling algo rithm 
MultiNest descr ibed by iFeroz fc HobsonI ((2QQ8[ ) and 
iFeroz et~al] (|2QQ9f ). We used 4000 live points with con¬ 
stant efficiency mode turned off, and set an enlargement 
factor of 0.1. The maximum a posteriori model param¬ 
eters and their associated 68.3% credible intervals are 
collected in Tables [ 8 ] and El Also found there are other 
derived parameters including the scaled semimajor axis 
(a/i?*), the orbital inclination angle i, the transit du¬ 
rations Ti 4 and T 23 (first-to-fourth and second-to-third 
contacts), and the planetary radius and semimajor axis, 
Rp and Op. Additional derived quantities of interest are 
discussed below. The folded transit light curves are dis¬ 
played graphically in Figure [TT] along with the models. 

The planetary radii for our targets are all under 2.4 


with the exception of KOI-2529.02, and the semimajor 
axes range from about 0.15 to 0.64 AU. We find that in 
all cases our planetary radii are larger than those cur¬ 
rently held on the NASA Exoplanet Archive. This may 
be partly due to the fact that our stellar radii reported 
earlier also tend to be larger, and that in four cases we 
have included dilution effects from close companions. 


[REST OE PAGE IS BLANK] 
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Table 8 

Final parameter estimates for the objects studied in this work. 


Parameter 

KOI-0571.05 

KOI-1422.04 

KOI-1422.05 

KOI-2529.02 

KOI-3255.01 

KOI-3284.01 

Fitted parameters 

P [days]. 

129.944y°;“« 

63.336001“;°°°50 

34.14204l°;°°°25 

64.00205l°;°°°g 

66.650621°;°°°“ 

op: oqqi q +0.00025 

on.zooiy_o 00029 

T [BJDutc - 2,400,000] .. 

55789.494o1o;oo4i 

55692.305i1“;°°26 

55686.01541“;°°“ 

55477.63071°;°°^ 

55670.68761°;°°“ 

55673.0307l°;°°44 

{Rp/R*) . 

0.0205inoi3 

0.0362tnM8 

u.uzy / _o 0037 

U.UC)O^_0 0035 

0.02791^042 

0 . 02001 °;°°“ 

log(p*,phot [gcm-^])"" . 

0.80lg;t? 

1.08l“j| 

0 77+0.44 

O'- ‘ '-0.19 

n 170 + 0.36 

O'-' ^- 0.21 

U.OU_o 19 

n 70+0.38 
‘^-0.28 

b . 

n 1O+0.28 

'^•^'^-0.13 

u.uy_o 09 

n 1 q+0.30 
*^•-*-^-0.18 

u.u^-o 09 

U.OD_o 31 

' 0 -^' 0 _o 34 

qi . 

U.C)O _0 28 

U.UO_o 03 

u.yo_o 18 

' 0 -^i-o. 2 i 

0 30+0-26 
U.aO-o 36 

D «3+0+5 

U.50-0 30 

q2 . 

U.UO_o 08 

0 03+0-22 
U.UO_o 03 

0 97+0-03 

0'-^' -0.15 

0 03+0-00 

u.ua_o 05 

u.oa-o 30 

n 34+0-32 

log/3 . 


(-0.532 ± 0.089)b 

(-0.532 ± 0.089)b 

(-1.330 ±0.084)l> 

(-0.080 ± 0.089)b 

(- 0.577 ± 0.070)‘> 


Other transit parameters 


{a/R^) 


ui . 

U2 . 

Ti 4 [hours] 
T 23 [hours] 


[^©] . 

ap[AU] . 

log(p^ [gcm-^])"" 

6min . 

'S'eff ['S'e] . 

'PpZ] [%] . 

7 ^ [rocky] [%] .... 


1781“ 

1371“ 

71 +29 
'1-10 

1041“ 

91+1 

7q + 26 

''^-15 

89.961°;“ 

89.95l°;“ 

SQ SQ + 0-11 

oy-oy_0.26 

89.93l°;“ 

89.901°;“ 

89.861°;“ 

0 70+0-69 

O'-' ’9_o.38 

0.35l°;“ 

1 p:9+0.32 

l-O^-0.34 

0-67l°;“ 

0'-^^-0.31 

0.88l°;tl 

—0 10+0-61 
'9-1’9-0.30 

— 0 Ol“*”0*-35 

0'-0'l_o.20 

—0 70+0-23 

O'- ''9_o.27 

— 0 14+0-38 
0'-l^-0.31 

— 0 31+0-25 
'^-'^^-0.33 

-o.42l°;“ 

0-’9^-0.24 

C).DU_q ig 

q f?9+0.26 

'l-0'^-0.26 

^•0"-0.34 

3 30+0-21 
0.09_q 29 

q 30+0.30 

'^•^^-0.41 

3 99+0.25 
^•^^-0.30 

3 26“^®'^® 
'^•^0)_Q 19 

3 34+0.31 

O.o^-O 26 

4 10+0-33 
"1-10-0.43 

3 91 +0-27 
^-^1-0.24 

3 24+0-61 

'^-^^-0.39 

Physical parameters 

1 1 7+0.08 
-0.08 

1 73+0-12 

1- '0>_o.l9 

1 45^ + 0-16 

1-40-0.25 

7q + 0.23 
''^-0.24 

2.14l°;“ 

1 19 + 0.16 
1-1^-0.17 

0 432+0.171 

O'•^'^^-0.053 

0.283l°;“° 

0 149+0-035 

0 33Q+0-164 
U.OOy_o 053 

0.2881°;°“ 

U.1DD_0 042 

0 734+0-022 

U. /O^-O 020 

0-915l°;°°° 

0 Ql 3+0.069 
u.yio_Q Qg4 

n 409+0.036 
*^-^^^-0.036 

0 504+0-035 
U.OU^_Q Q35 

n 71 7+0.079 
O'-'1'-0.156 

p 

b 

1 + 
0 0 
b b 

0 12+0-10 
'^•-'■^-0.09 

O'-IO'-O.IO 

0 10+0-13 
O'-l^-O.ll 

0'-’4^-0.02 

U.U9_q Q3 

0.30l°;l° 

0 39+0-10 

u.o^-o 16 

1 45+0-60 

i.oy-o 79 

9 1 3+0-74 
^•1^-0.88 

1 40+0-67 

1-^'0_Q 77 

98.4 

99.7 

82.0 

72.3 

56.8 

71.8 

68.4 

30.6 

50.7 

0.6 

11.7 

69.6 


Note. — All quantities correspond to the mode of the posterior distributions. 
^ Mean stellar density derived from the light curve fit. 

^ Gaussian prior. 

^ Mean stellar density derived from our stellar evolution modeling. 
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Table 9 

Final parameter estimates for the objects studied in this work. 


Parameter 


KOI-4005.01 KOI-4087.01 KOI-4427.01 KOI-4622.01 KOI-4742.01 KOI-4745.01 


P [days]. 

T [BJDutc - 2,400,000] .. 

(Rp/R*) . 

log(/9*,phot [gcm“3])a . 

b . 

91 . 

92 . 

log/3 . 


Fitted parameters 


178.13961°;“!® 

101 1 1 141 +0-00087 

iUi.iii^i_0 00068 

147.66061°;°°!! 

55399.3987l°;oog4 

55756.2986l°;oo2? 

55815.26161°;°°!^ 

'^•’^^'^^^-0.00111 

0.03038l°;““2 

0.0339in!l2o 

0 99+0-16 

0.251°;!° 

0.78l°;°! 

'^•-‘-'^-0.10 

u.uy_o 09 

0-08l°;^® 

U.UO_o 05 

0 03+0-10 
U.U9_o 03 


U.UO_o 08 

0 03+0-20 
U.U9_o 03 

U.UO_o 08 


(-0.804 ± 0.089)1> 


207.2482l°;°“2 

112 30'53+0-0024 

11Z.OUOO_0 0028 

177 

1' /•ODy9_o 0030 

55667.77641°;°°!° 

55849.55781!'°“! 

55630.24601!;!!“ 

0.0280ino!I 

0.02111!;“!! 

u.uou9:_o 0022 

0.89l°;“ 

0.68l°;“ 

n qo+0.30 
'^•'^^-0.20 

0.25l°;°° 


'0-^O_0 26 

0 03+0-15 
U.U9_o 03 

0 99+®-®® 

n qq + 0.32 

U.OO_o 33 

U.U9_o 03 

u.yu_o 30 

0 1 9+0-29 


(a/i?*) 


ui . 

U2 . 

Ti 4 [hours] 
T 23 [hours] 


[^©] . 

ap[AU] . 

log(P3^ [gcm-^])"" 

6min . 

'S'eff ['S'e] . 

V[RZ] [%] . 

7 ^ [rocky] [%] .... 


Other transit parameters 


1421!® 

99.o1!°42 

190l“ 

260l“ 

1461“ 

i5il!! 

89.951!;“ 

89.931!;“ 

89 97+®-®® 

ov.iyi _0.08 

89.971!;“ 

89.941!;“ 

89.941!;“ 

0.391!;!! 

0.341!;!! 

0 07“^®'^^ 

-0.39 

0.281!;“ 

1 03”^®'®® 

i.uo_o 49 

0 65+®-^^ 

U.DO_o 48 

-0.081!;“ 

-0.051!;“ 

—0 23“*”®-^® 
'^•^'^-0.34 

—0 02“*”®'®® 
'®-'®^-0.20 

—0 43+®-®2 

-0.151!;!! 

^•"^^-0.42 

/ •y9_o 21 

5 99+®-®® 
"^•^^-0.29 

6 05”*”®'®® 

D.UD_o 32 

5 62+®'^® 
®-’^^-0.52 

S'SSI!;!! 

^•*0^-0.35 

<-7 r\ij-\-0.23 
-0.33 

5 42“*”®-^® 
0-4^_0.37 

5 01”*”®'®® 
^^-0.32 

5 26”*”®'^® 
®-^'^-0.52 

0 (^o+0.70 

®-^®-0.56 

Physical parameters 

2 24+0-16 
^•^^-0.45 

-’-•^^-0.19 

1 84+®-^^ 
-*-•^^-0.24 

1 04+®-22 
-'■•'^^-0.24 

1 34+0.11 
-'■•'^^-0.18 

9 qq+0.19 

z.oo_o_22 

U.ODC)_o 080 

0 242“^®-®®® 

U.Z4Z_o 041 

0 419+®-®^® 
u.^±j^_0 081 

0 64”^®'®^ 

0 409”^®'^®® 

U.^U»_Q QgQ 

0 495^®'^®® 

u.ayo_Q 

u.zyo_o 076 

0 o37+®-®®4 

U.DO/ _o 162 

n +0.083 

U. /Oi_o 095 

0 672+®-®®^ 
U.D/Z_o 211 

0 005”*”®'®®® 
U.DUO_o 035 

0 477+®-®43 
' -0.035 

0 03+®-®® 
U.U9_o 03 

0 34+0-12 

’®-’®^-0.02 

0 10”^®'^^ 
O'-^’O-O.IO 

0 04”*”®'®® 

'®-'®^-0.04 

O-lllo.n 

1 S3+®-51 
l.OO_0 62 

1 20“*”®'^® 
^•^'0-0.65 

n 000+0.069 

'4-^oO_q iiQ 

0 21”*”®'^^ 

0 66”*”®'^® 
U-OD_o 41 

0.861!;“ 

74.1 

89.0 

71.9 

48.6 

96.9 

89.9 

6.5 

29.8 

27.3 

45.0 

60.7 

4.9 


Note. — All quantities correspond to the mode of the posterior distributions. 
^ Mean stellar density derived from the light curve fit. 

^ Gaussian prior. 

^ Mean stellar density derived from our stellar evolution modeling. 
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8. CONSTRAINTS FROM ASTERODENSITY PROFILING 

A comparison between the mean stellar density in¬ 
ferred from our light-curve fits above (photometric den¬ 
sity p^,phot) and some independent measure of the den¬ 
sity, such as that derived from model isochrones, can pro¬ 
vide very useful information about a planetary system. 
This technique, referred to as Asterodensity Profil i ng 
(AP), was discussed in detail bv lKipping et al.l (|2Q12[) Fq 
For a “vanilla” exoplanet (circular Keplerian orbit, no 
blends, etc.) these two density metrics should agree. De¬ 
viations between them occur when one of the following 
physical effects is not accounted for in the light curve 
analysis: i) orbital eccentricity (‘photo-eccentric effect’); 
ii) contaminating light (‘photo-blend effect’); in) unoc¬ 
culted star spots (‘photo-spot effect’); iv) transit tim¬ 
ing variations (‘photo-timing effect’); v) transit duration 
variations (‘photo-duration effect’); vi) substantial mass 
of the transiting object (‘photo-mass effect’); or vii) the 
object orbits a differe nt star. These e ffects have been 
described in detail by iKippind (|2Q14a[) . who provided 
approximate expressions for each. 

In what follows we apply AP to 1) test the hypothesis 
that all planets in the multiple transiting planet systems 
KOI-0571 and KOI-1422 orbit the s ame star (Multi-bod y 
Asterodensity Profiling, or MAP; iKipping et al.|[2QT^ : 
2) derive minimum eccentricities for the twelve KOIs in 
our sample; and 3) evaluate the possibility that the val¬ 
idated planet orbits the close companion instead of the 
target, in the four cases where high-resolution imaging 
has identified such sources. 

8.1. Using AP to derive average stellar densities for the 
multi-planet systems KOI-0571 and KOI -1422 

Of the twelve candidates studied in this work, three 
reside in systems with five transiting planets each: one 
orbits KOI-0571, and two are associated with KOI-1422. 
These rare systems provide extra information that we 
will exploit here. For each planetary candidate in each 
system we performed an independent detrending and 
subsequent light curve fit using the methods described 
earlier. 

If all five planets in each of these two systems orbit the 
same star, then we expect all of the planets to maintain 
low orbital eccentricities in order for these compact sys- 
tem s to be dynamica lly stable over Gyr timescales (see 
also lRowe et al.l[2 01 41) . Therefore, under the assumption 
that the planets orbit the same star, one would expect 
negligible photo-eccentric effects, implying that the mean 
stellar densities derived from the light curves should be 
consistent. On the same assumption that all five planets 
orbit the same star, two other AP effects listed above 
can also be neglected: the photo-blend effect, and effects 
from the planets orbiting different stars. The latter sce¬ 
nario is obviously not relevant under the posed assump¬ 
tion, and the former can be negated since, to first order, 
contaminating light that is unaccounted for disturbs the 
light-curve-derived stellar density for all five planets to 
the same degree. 

A comparison of the photometric stellar densities for 
the five planets around KOI-0571 and KOI-1422 reveals 

A similar concept was described bv lTinglev et all (|20Iin . al¬ 
though the current implementation is considerably more advanced. 


excellent agreement in both cases, as shown in Figure [121 
Furthermore, similar diagrams (not shown) comparing 
the limb-darkening parameters Ui and U 2 also show very 
good agreement among all planets in both systems. We 
conclude from this that each system is consistent with 
having five planets orbiting the same star, given the avail¬ 
able data. 

Averaging over the four other planets in KOI-0571, 
and the three others in KOI-1422, we may derive a 
mean photometric density for the host stars that we de¬ 
note p^,map- We obtained p*,map = 5.4llQ;24gcm“^ 
and ^ respectively, where the extra dilution 

from the close companion to KOI-1422 has been properly 
accounted for, as before. The la and 2cr confidence inter¬ 
vals of these values are marked in Figure [12] with dashed 
lines. Both of these mean densities are consistent with 
the corresponding photometric densities derived from 
the light curves of KOI-0571.05, KOI-1422.04, and KOI- 
1422.05 individually, which are p*,phot = 5.9I2.5 gcm“^, 

ll.dl^^Q gcm“^, and 5.3l2!4 gcm“^, respectively. 

The average stellar densities for these multi-planet sys¬ 
tems (computed from the remaining planets in each case, 
excluding the validated planets) were used as luminosity 
indicators to strengthen the determination of the host 
star par ameters via our ste llar evolution modeling in 
Sect. [5] (|Sozzetti et al.l [20071) . This is particularly use¬ 
ful for these three KOIs as neither of the host stars has 
a spectroscopic determination of log g available to other¬ 
wise constrain the luminosity (see Table [5]). 

8.2. Using AP to measure minimum eeeentrieities 

A measure of the stellar density that is independent 
of that derived from the individual light curves of the 
validated planets is available for each of the twelve 
KOIs studied in this work (see Sect. [5]). For KOI-0571 
and KOI-1422, this comes from the MAP-based density 
(p*,MAp) refined by our isochrone analysis by making 
use of additional constraints on age from gyrochronology, 
along with the spectroscopic temperature and metallic- 
ity estimates. For KOI-3284, 4087, 4427, and 4622 our 
stellar evolution modeling used only the age along with 
Teff and [Fe/H], and for the remaining targets that are 
earlier in spectral type we used log^ as the luminosity 
indicator, which for these stars we were able to determine 
spectroscopically using SpecMatch. 

These independent stellar densities, denoted here sim¬ 
ply by may be compared with the light-curve-derived 
stellar density of each validated planet, p*,phot- In these 
cases we expect the photo-eccentric effect to be the dom¬ 
inant AP effect since contamination has been accounted 
for in the fits (photo-blend effect), the stars show rel¬ 
atively low activity (photo-spot effect), there are no 
known timing effects for these objects (photo-timing and 
photo-duration effects), and their small sizes suggest they 
should be of low mass (photo-mass effect). We therefore 
attribute any differences between p*,phot and as being 
due to orbital eccentricity. The minimum eccentricity 
for each KOI is directly given by t his com parison via the 
expression presented bv iKippina (j2Q14a[ ). 

„ _ |l-(p*,phot/p*)^/"l 

l + (p*,phot/p*)^/^ ■ 
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Figure 11. Light curve fits and i^ep/er observations (small dots) for our twelve KOIs. Pink circles represent 10-point binned data. 


The results of this calculation, shown in Tables [8] and [9l 
reveal that most of the KOIs in the sample are consistent 
with orbiting their parent stars on low-eccentricity orbits. 

8.3. Using AP to test blend seenarios 

For the four KOIs found to have close companions from 
our high-resolution imaging (KOI-1422, 2529, 3255, and 
3284) the same techniques from the preceding section 
enable us to revisit the possibility that the planets orbit 
the companions rather than the primary stars. As done 
earlier, we made the reasonable assumption that these 
close neighbors are bound to the targets. This allows us 
to reuse the target isochrones to infer approximate stel¬ 
lar properties for the companions, including their mean 
densities, on the basis of their brightness differences com¬ 
pared to the target. We re-fitted the transits adopting re¬ 
vised blend factors (/3) appropriate for the different host 
stars to derive photometric densities for the companions, 
P*,phot- We then used these densities to estimate a min¬ 
imum eccentricity for the scenarios in which the planets 


orbit the companions rather than the brighter stars. 

We report the results of this calculation in Table [TOl 
where we list also the planetary radii that we infer if 
the planets transited the companions, as well as their ef¬ 
fective insolation level (see next section). Although the 
eccentricity distribution of s mall planets arou nd M dwarf 
stars is presently unknown (jKippindl2Q14b[ ). in general 
we expect circular orbits to be more likely than eccen¬ 
tric orbits, based on evidence from ot her sub-populations 
of planets (see, e.g., lKippingll2Q13al) . We may thus use 
emin as a rough criterion to decide which location for 
the planet is favored. Additionally, the sign of the quan¬ 
tity log(p*^phot/p*) conveys useful information. From ge¬ 
ometric arguments, a planet in an eccentric orbit has 
a higher prob ability of tra nsiting near periapsis than 
apoapsis (e.g.. lBarnesll2QQ7l ). For the more probable case 
of 0 < cj < TT (“near” periapsis transits), one expects 
log(p^^phot/p*) > 0- la contrast, for the less probable 
scenario of tt < cj < 27r (“near ” apoapsis trans its), one 
expects log(p*^phot/p*) < 0 (see lKippingll2Q14al ). There- 
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Table 10 

AP-derived minimum eccentricities and other properties for 
KOIs with close companions. 


Parameter 

Primary 

Companion 

KOI-1422.04 

Kp [mag]. 

16.15 ±0.03 

17.72 ±0.15 

log/3. 

(-0.628 ± 0.061) 

(0.628 ± 0.061) 

log(p^ [gcm-3]). 

0 .yiO_Q Qg 4 

l-342l°;0|^ 

log(p*,phot [gcm-3]) .. 


l-O+ols 

Iog(p*,phot/P*). 

0 .iD_Q 

-0.14 

^ • 

0 . 12 +°-nQ 

0 . 20 +n on 



— 0.20 

[^©1 . 

1.76+°;}8 

^•^1-0.18 

‘S'eff [S^] . 

0.38+°“ 

n 90 + 0 .O 6 
*^•^^- 0.10 

KOI-1422.05 

Kp [mag]. 

16.15 ±0.03 

17.72 ±0.15 

log/3. 

(-0.628 ± 0.061) 

(0.628 ± 0.061) 

log(p^ [gcm“3]). 

u.yio_Q Qg 4 

1 . 0 UZ_Q Q 3 g 

log(p*,phot [gcm-3]) .. 

0.^^_Q4Lg 

0.77l°;|| 

Iog(p*,phot/P*). 

—0 14+0-^® 

-0.56l°;« 

p • 

0 . 10 +nin 

0.79+n 


— 0.10 

^ — 0.22 

Rp [-R©]. 

1 42+0-19 

1 . 68-0 53 

‘S'eff [S^] . 

qc:+0.49 
1 . 00-0 66 

1 ni+0-50 
J--'01-o.47 

KOI-2529.02 

Kv [mag]. 

15.86 ±0.10 

22.00 ± 0.20 

log /3. 

(-1.330 ±0.084) 

(2.474 ± 0.089) 

log(p^ [gcm-3]). 

n /loo+o.oso 

o.48^_q Q 3 g 

1-8661^39 

log(p*,phot [gcm-^]) .. 

n ^To + O.SO 

O./Z_Q 2 i 

0 . 88 l°;“ 

Iog(p*,phot/P*). 

0 24+®-'^® 
U.ZU_Q 21 

—1 01 +°+® 

^ • 

0 iq+0-13 

0-65l°;°° 


Rp [^©]. 

P f^q + 0.23 
^•''^-0.24 

8 - 21^1 

‘S'eff [‘S'e]. 

i.oy-o 79 

0.3o1°;°8 

KOI-3255.01 

Kp [mag]. 

15.01 ±0.10 

15.21 ±0.20 

log/3. 

(-0.080 ± 0.089) 

(0.080 ± 0.089) 

log(p* [gcm-3]). 

U.OU4:_q Q 3 g 

U. 0 ^ 4-0 010 

log(p*,phot [gcm“^]) .. 

U.OU_Qj^g 

0 51+0-16 

'0-^l-0.18 

Iog(p*,phot/P*). 

0 01 + 0-22 

U.UU -0 16 

p . 

0.02 + n n? 

0 . 02 +n n? 


— 0.02 

— 0.02 

Rp [-^©]. 

2 14+®-^^ 
^•1^-0.17 

9 9 [;:+ 0.20 
^•^"^- 0.20 

8 'eff [‘S'e]. 

9 1 c;+0-r4 
^• 1 ^- 0.88 

1 Q1+0-57 

J--^J--0.49 

KOI-3284.01 

Kv [mag]. 

14.58 ±0.10 

17.01 ± 0.20 

log d. 

(-0.577 ± 0.070) 

(1.036 ± 0.088) 

log(p* [gcm-3]). 

n 71 7+0-0’^9 

'0./i/_Q4L5g 

-| 091 © 0.020 

i.9Zi-o oil 

log(p*,phot [gcm-^]) .. 

n 70 + 0-38 
^•'^-0.28 

0.77l°;« 

Iog(p*,phot/P*). 

U.UO_Q 35 

-0.54l°;« 

p ■ 

0.03 + n nS 

0.40+no? 


— 0.03 

— 0.31 

Rp [^©1. 

1 19 + 0.16 

1 77+0.23 

1 - ‘ ‘ -0.21 

'S'eff [Sq] . 

1 40+0-®’^ 

0 -86i°i° 


Note. — The log (3 values indicated in parentheses were used as 
Gaussian priors. For KOI-2529.02 and KOI-3284.01 the /3 factors 
for the companions are not the reciprocal of the values for the 
primaries because of the presence of additional companions in the 
aperture (see Table [2]). Those additional stars are ruled out a s the 
source of the transit signals by our centroid analysis of Sect. 14.21 
but still cause dilution. The symbol p* represents the mean stellar 
density based on our isochrone analysis (and other constraints), 
while p*,phot is the photometric density from the individual light 
curve fit for the planet. 




Figure 12. Mean stellar densities from separate light curve fits 
to each of the planets in the five-planet systems KOI-0571 and 
1422. The dashed lines and gray areas represent the Icr and 2cr 
confidence regions of the average density computed from the other 
companions in each system (open circles) excluding the ones we 
have validated here (filled circles). In both cases the agreement 
between this mean density and the separately determined one of 
the validated planets is excellent, supporting the notion that in 
each system all five planets orbit the same star. 

fore, in general the photo-eccentric effect is more likely 
to (but will not necessarily) produce a positive sign for 
this logarithmic ratio. Additionally, a diluted transit 
in a circular orbit will cause lo g(p^.nhot/p^) to be nega¬ 
tive via the photo-blend effect (jKipping|l2Q14aI ). but can 
never yield a positive sign. For these reasons, finding 
log(p^^phot/p*) > 0 is more consistent with a bona-fide 
planet than a negative sign. 

As seen in Table [TOl despite the large uncertainties the 
minimum eccentricities are considerably higher if KOI- 
1422.05, KOI-2529.02 and KOI-3284.01 orbit the sec¬ 
ondary star, which seems a priori less likely, although 
quantifying an accurate odds ratio is not possible without 
knowing the eccentricity distribution of terrestrial-sized 
planets around M dwarfs. For the first two cases this 
assessment agrees with our conclusion from the BLENDER 
study (Table [7j), which also favored the planet transiting 
the primary. For KOI-2529.02 and KOI-3284.01 the sign 
of the log-density ratio log(p*^phot/p*) is also negative 
if the planet is around the companion, requiring either 
geometrically disfavored near-apoapsis transiting plan¬ 
ets, or simply the photo-blend effect. For KOI-1422.04 
the minimum eccentricities between the two scenarios are 
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comparable, but the sign of the log-density ratio favors 
the primary star as the host. Finally, the test for KOI- 
3255.01 is ambiguous according to both metrics, in agree¬ 
ment with the indications from BLENDER that were also 
inconclusive for this planet. 

9. HABITABILITY 

While a broadband transit light curve alone does not 
reveal the “habitability” of a given planet, it does allow 
us to determine how close to the “habitable zone” a par¬ 
ticular planet is. Therefore, an important distinction to 
make is that we do not claim any of the planets in our 
sample are truly habitable, but merely that the insola¬ 
tion they receive from the host star is suitable for the 
presence of liquid water on the surface across the range 
o f plausible atmospheric c onditions. 

iKopparapn et ahl (|2Q13[ ) have presented calculations of 
the habitable zone for a wide variety of atmospheric con¬ 
ditions, and considered also empirical boundaries for the 
hot inner edge modeled by a recent-Venus-like world, and 
a cold outer e dge modeled by a n early-Mars-like world. 
More recently IZsom et ahl (j2Q13[ ) took a broader view of 
the conditions that might provide habitable surface con¬ 
ditions closer to the host star, and postulated an inner 
boundary for the HZ that can be as close as 0.38 AU for 
a solar-type star if the relative humidity of the plane¬ 
tary atmosphere is low and the albedo high. Given the 
complexity of the problem and our still limited state of 
knowledge, in this work we have chosen to adopt the 
broadest possible limits that can plausibly lead to suit¬ 
able conditions for life, in order to avoid prematurely 
dismissing as uninteresting planets that may yet be hab¬ 
itable. Specifically, we adopted the empirical c old outer 
boundary for the HZ of iKopparapn et ahl (j2Q13[ ) modeled 
by an early-Mars- like world, and the inner edge for dry 
desert worlds from IZsom et al.l (j2Q13f ). These boundaries 
depend upon two key input parameters: S'eff, the effec¬ 
tive insolation received by the planet, and the spectral 
type of the star, for which the usual proxy is Teff, the 
star’s effective temperature. We define 5'eff in terms of 
the Earth’s insolation, S'0, as 

^eff ^'kj Lq 

where is the stellar luminosity and the planetary 
semimajor axis in astronomical units. For each planet 
in our sample we computed a joint posterior for S'eff and 
Teff based on the stellar properties inferred in Sect.[5l as¬ 
suming the planet’s orbit is circular. At each realization 
in this posterior we evaluated whether the sample falls 
inside or outside the adopted HZ boundaries. Count¬ 
ing the total number of cases inside allows us to assign 
a statistical confidence that the planet in question is in 
the HZ, denoted T[HZ]. We provide those results in Ta¬ 
bles [8] and [9l This statistic al approach to the h abitable 
zone was first described by iKipping et al.l (j2Q13f ) for the 
planet Kepler-22 b. 

Additionally, we consider the radius of the planet 
to be an important parameter in assessing its poten¬ 
tial for habitability, as this can determine whether 
the planet is likely to have a solid surface that 
can support liquid water. In the past it was of¬ 
ten conventionally assumed that there likely existed 



Figure 13. Joint posterior distribution of the planetary radius 
(y axis) and effective insolation {x axis) for KOI-057I.05 (Kepler- 
186 f). Labeled vertical lines ma rk the various boundari es of 
the habitable zon e as defined by IKopparapu et aP (120131) and 
IZsom et ahl (120131) . Here we adopt the optimistic ‘dry desert’ and 
’early Mars’ inner and outer edges, respectively. These lines are 
drawn assuming an effective stellar temperature equal to the mode 
of the associated posteriors. Contours enclose 68.3%, 95.4% and 
99.7% of the Monte Carlo samples (corresponding to the Icr, 2(7, 
and 3(7 levels). 

some boundary in radius between ro cky planets and 
mini-Nept nnes (e.g.. IWeiss et al.ll2Q13[) . perhaps around 
1 .75^0 (|LoDez fc FortncA^ I2Q13D . However, recent 
counter-examples to this paradigm make such a sim¬ 
plistic division questionable. For exam ple, KOI-314c 
is a rn ini-Neptune despite being 1.6 (IKipping et al.l 
l2Q14aD . and Kepler- 10c has been claimed t o be rocky de¬ 
spite bei ng 2.35 T0 (jPuninsaue et al.ll2Ql4 for a d iffering 
view see lRogersll2Ql3 and IWolf gang fc Lopezll2Q14D . Nev¬ 
ertheless, it is reasonable to assume that smaller planets 
are more likely to be rocky, and thus we present joint pos¬ 
terior distributions in radius and insolation. An example 
of this diagram for KOI-0571.05 is given in Figure [T3l in 
which we indicate various boun darie s for the HZ as de - 
fined bv iKopparapu et ^ (|2Q13[) and lZsoni et al.l (|2Q13[) . 
Similar graphs for the remaining targets in our sample 
are shown in Figure [TTl 

From the standpoint of habitability it is of interest 
also to quantify the likelihood that the planets are rocky, 
given their measured sizes and corresponding uncertain¬ 
ties. To do this we have considered an empirical den¬ 
sity model giving a statistical estimate of the mass of 
a planet as a function of its radius, constructed based 
on the available sample of transiting planets with mea¬ 
sured densities (i.e., with masses measured either via the 
Doppler technique, or from TTVs). This model is ex¬ 
pected to be biased towards larger and more massive 
planets that are easier to detect, but it is possible to re¬ 
move the bias by convolving the density model with es¬ 
timates of th e rate of occurrence of planets as a function 
of their size (|Fressin et al.ll2Q13D. which we have done in 
a Monte Carlo fashion (|Fressin et al.ll2Q15l) . In this way 
we may infer the true distribution of planets in the mass- 
radius diagram, which can then be used to establish the 
fraction of rocky or gaseous planets at a given radius or 
mass. For the purposes of this work we have considered 
any planet denser than a composition of 100 % MgSiOs 
perovskite to be rocky (|Seager et al.ll2QQ^ . Using the 
above model and our posterior distributions for the ra¬ 
dius of the planets in our sample, we have computed the 
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Figure 14. Similar to Fig. ll3l for the other planets in this study. 

probability that they are rocky, 7^[rocky], by assigning a 
mass to each link in the chain and counting the number 
of samples denser than the threshold, analogous to what 
was done previously for the HZ. The 7^[rocky] values are 
listed for each object in Tables [8] and [9l Estimates of 
7 ^ [rocky] using an altern ate density model proposed re¬ 
cently by iRogersI (|2Q14f ) lead to similar assessments, at 
least qualitatively. 

The results of our calculations indicate that all of the 
planets in our sample have a ^50% or greater chance 
of being in the temperate region of their parent stars. 
KOI-0571.05, 1422.04, and 4742.01 are the best candi¬ 
dates for residing within the HZ boundaries, with 7^[HZ] 
equal to 98.4%, 99.7%, and 96.9%, respectively. The first 
two receive roughly one third of the flux of the Earth, 
while the latter receives two thirds. Based on their mea¬ 
sured sizes, KOI-0571.05 and KOI-4742.01 have the best 
chance of being rocky; KOI-1422.04 is larger and has a 
somewhat smaller chance. KOI-3284.01 is the object in 
our sample with the highest probability of being rocky 


in nature (69.6%); its size is consistent with being the 
same as that of the Earth: Rp = R^. Other 

planets with non-negligible chances of bein g r ocky are 
KOI-1422.05, 4087.01, 4622.01, and 4742.01H On the 
other hand, the larger planets KOI-2529.02, 4005.01, and 
4745.01 have 7^[rocky] ^ 5% or smaller, and are most 
likely mini-Neptunes. KOI-3255.01 may still be rocky 
despite its radius of 2.14 R0, with 7^[rocky] = 11.7%. 

In terms of their insolation, four of our validated plan¬ 
ets with a sizable chance of being rocky have values 
within 50% of that of the Earth: KOI-1422.05, 3284.01, 
4087.01, and 4742.01. Of these worlds, KOI-3284.01 and 
KOI-4742.01 appear to be the most “Earth-like” when 
considering both their size and amount of incident flux 
jointly, surpassing KOI-0571.05, which is also rocky but 
receives only one third of the incident flux of our planet. 

KOI-4427.01 would otherwise join this group, but it is not 
formally validated as a planet to the same level as the others. Its 
insolation is only about a quarter of that of the Earth. 
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10. DISCUSSION AND CONCLUDING REMARKS 

Eleven of the twelve candidates in our sample have 
been validated here as planets to a very high degree of 
confidence exceeding 99.73% (3cr). The high bar we have 
set in this work is commensurate with the importance 
we assign to these planets, which are all small and po¬ 
tentially in the HZ of their host stars. Few such objects 
have been verified so far, and our work now essentially 
doubles the number of cases listed in Sect. [H 

Three of our targets, KOI-0571.05, KOI-1422.04, and 
KOI-1422.05 are in systems of five planets each, and have 
been previously validated by others using very different 
methodologies and assumptions. Our BLENDER valida¬ 
tions confirm those assessments using conservative and 
realistic hypotheses, generally with a higher degree of 
confidence. KOI- 1422.04 and KOI-1422 .05 w ere included 
in the studies of iLissaner et al.l (|2014[) and iRowe et ahl 
(|2014[ ). who considered them to be validated based on 
the general conclusion that the vast majority (> 99%) 
of Kepler candidates in multiple systems are lik ely to 
be bona-fide planets, as they showed statistically^ In 
the present work we are able to assign a specific con¬ 
fidence level (99.74% and 99.89%) to these two plan- 
ets using BLENDER. KOI-Q 571.05 was validated separately 
by I Quintana et al.l (|2014[) through statistical arguments 
somewhat similar to those in BL ENDER in some re s pects , 
but not in detail. For example, I Quint ana et al.l (|2014f ) 
included short period planets in their estimate of the 
planet prior, whereas we chose to count only planets with 
similar period as the candidate, because the rates of oc¬ 
currence of transiting planets are a strong function of 
orbital period. For a planet with a long period of ^130 
days such as KOI-0571.05, including short periods over¬ 
estimates PL by about a factor of 50 according to our 
calculations. On the other hand, including short periods 
als o for the false positiv es (both for planets and EBs), 
as iQnintana et al.[ (|2Q14D have done, will tend to make 
the blend frequencies higher. This partly offsets the bias 
from the larger planet prior, but still results in a final 
odds ratio a factor of 2 greater than if short periods are 
not included. An additional difference is that in order 
to rea ch a comfortable validation level IQnintana et al.l 
(|2Q14[ ) applied a “multiplicity boost” to their planet prior 
(increasing it by a factor of 30), invoking the fact that 
KOI-0571.05 is in a five-planet system, that the other 
four p lanet s were considered v alidated by ILissaner et al.l 
(|2Q14[ ) and iRowe et al.1 (|2Q14D . and that false positives 
are significantly less common in mult is. While this is 
a perfectly legitimate argument, the odds ratio achieved 
with BLENDER is high enough (for all our multis) that this 
boost was not required here. Had we used it, our confi¬ 
dence level for KOI-0571.05 would be 30 times higher. 

One of our candidates, KOI-4427.01, does not quite 
reach our threshold for validation, although in other 
contexts a 99.16% confidence level such as we achieve 
might be considered high enough to declare it a bona- 
fide planet. Interestingly, this is the same candidate for 
which our centroid motion analysis indicated less confi¬ 
dence than in other cases that the transit signal comes 

We note that another of our targets, KOI-2529.02, is a candi¬ 
date in a two-planet system that was also included in the statistical 
studies above, but it was not considered to be validated there due 
to unresolved issues with the follow-up observations. 


from the target itself as opposed to a nearby background 
location, consistent with indications from BLENDER. We 
note also that the host star does not show a clear ro¬ 
tational signature, which causes its properties to be less 
well determined than others for which we were able to 
use an age estimate from gyrochronology to strengthen 
the model fits. For these reasons we do not include this 
object among our formally validated planets, although it 
is entirely possible that additional follow-up observations 
may enable a higher level of validation in the future. It 
remains a very interesting object, as it is likely to be 
within the HZ and also has a fair chance of being rocky. 

We find that the mean densities derived photomet¬ 
rically from the individual planets in each of the five- 
planet systems KOI-0571 and KOI-1422 show very good 
agreement, supporting the notion (for the first time in 
KOI-1422) that in each case the five planets orbit the 
same star. This is further supported by a similar con¬ 
sistency found among the limb-darkening coefficients {ui 
and 1^2), which were fitted for separately for each of the 
planets. For the four targets in our sample with close 
companions the location of the planets is more uncer¬ 
tain, however. In one case (KOI-2529.02) both BLENDER 
and AP strongly favor the planet transiting the primary 
star. For KOI-1422.04, KOI-1422.05, and KOI-3284.01 
the evidence is in the same direction but not as strong, 
and for KOI-3255.01 both BLENDER and AP are inconclu¬ 
sive as to which star has the planet. In these four cases 
even if the planet orbits the companion the inferred radii 
are still relatively small (1.68-2.25 770), and at least for 
KOI-1422.05 and KOI-3284.01 there is still a good chance 
they are rocky. 

Except for KOI-2529.02, the inferred radii of all our 
validated planets are nominally smaller than 2AR^. In 
some cases they are rather different than those reported 
previously, either because of differences in the stellar 
radii, or because some of the previous studies have not 
accounted for dilution by nearby companions. The stellar 
properties in this work are all based on spectroscopic esti¬ 
mates of the temperature and metallicity, as well as log g 
for the hotter stars and other constraints on luminos¬ 
ity for the cooler ones, superseding properties reported 
previously that have relied only on photometry. 

This study has significantly increased the number of 
small validated transiting planets (< 2.5 770) that are 
potentially in the HZ of their host star. Excluding 
those in our sample that have been announced pre¬ 
viously, as well as KOI-4427.01, we now add to the 
list KOI-1422.05(*), KOI-2529.02, KOI-3255.01(*), KOI- 
3284.01(*), KOI-4005.01, KOI-4087.01(*), KOI-4622(*), 
KOI-4742.01(*), and KOI-4745.01. The ones followed by 
an asterisk are small enough to be rocky, with a some¬ 
what reduced chance of that being the case for KOI- 
3255.01. We note that rather than relying on an esti¬ 
mate of the equilibrium temperature for our assessments 
regarding the HZ, here we hav e used the effective insola¬ 
tion, which, as pointed out bv iKopparapn et ahl (|2Q13D . 
dispenses with the uncertainty in the albedo. Our HZ 
statements are quantified by providing a probability that 
the planet lies within the adopted inner and outer edges 
of the region, which properly accounts for the uncer¬ 
tainties in all measured stellar and planetary properties. 
Similarly, we have provided a quantitative indication of 
the probability that each of these planets is rocky, also 
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accounting for uncertainties. 

A number of the previously announced small HZ plan¬ 
ets (see Sect. [T]) have radii in excess of 2 i?0, and accord¬ 
ing to our prescription for 7^ [rocky] they are significantly 
less likely to be rocky than perhaps deemed to be the case 
in the original publications (particularly Kepler-22 b and 
Kepler-296f, with radii of 2.38i?0 and 2.31 i?0, respec¬ 
tively). Thus, with the new examples added here, the list 
of rocky planets in the HZ is enlarged by a factor of two. 
KOI-3284.01 and KOI-4742.01 are now the validated, 
rocky, HZ planets that appear most similar to the Earth 
when considering both their size (Rp = and 

Rp = 1.34lo!i8^©) and insolation (5'eff = l-4Olo.77 5'0 
and Seff = 0.66to;4f ^e) jointly. 
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